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PREFACE 

This  report  contains  descriptions  and  operating  instructions 
for  a  collection  of  general  utility  programs  applicable  to  various 
phases  of  crystal  structure  analysis  with  computers.  All  coding 
is  done  in  FORTRAN  II  language  for  32K  core  IBM  704  and  7090 
machines;  the  FORTRAN-language  programs  are  completely  listed 
in  the  Appendix. 

The  data  input  to  many  of  the  programs  has  been  designed  to 
be  compatible  with  the  output  of  the  automatic  data  reduction  pro¬ 
gram  XRDDR  (H.G.  Norment,  “An  X-Ray  Diffraction  Data  Reduc¬ 
tion  Program  for  the  IBM  704  and  7090,”  NRL  Report  5739,  Feb. 
1962). 

Where  data  input  from  special  tape  is  required,  the  logical 
tape  numbers  are  assigned  symbolic  designations  near  the  begin¬ 
nings  of  the  FORTRAN  programs.  Thus,  if  it  is  required  that  these 
tape  numbers  be  changed,  it  is  sufficient  to  change  only  the  one 
FORTRAN  statement  card  and  then  recompile. 

In  the  form  presented  in  this  report,  the  programs  are  writ¬ 
ten  for  use  with  the  IBM  7090  IB  Monitor  system. 


PROBLEM  STATUS 

This  is  a  final  report  on  one  phase  of  the  problem;  work  on 
other  phases  continues. 


AUTHORIZATION 

NRL  Problem  C07-03 
Project  RR  001-02-43-4805 


Manuscript  submitted  November  21,  1962. 
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A  COLLECTION  OF  FORTRAN  PROGRAMS 
FOR  CRYSTAL  STRUCTURE  ANALYSIS 


1.  UNIT  CELL  PARAMETERS  AND  ERRORS  BY  LEAST  SQUARES 

1.1  Purpose 

Given  a  set  of  reflection  data  consisting  of  Miller  indices  and  sin  B  or  sin  e/\  values, 
compute  the  least-squares  reciprocal  and  real-space  unit  cell  parameters  and  uncertain¬ 
ties  for  any  crystal  system. 

1.2  Input 

The  data  input  is  taken  from  cards  immediately  following  the  program  deck.  It  con¬ 
sists  of  one  control  card  and  a  deck  of  reflection  cards. 

A.  Control  Card: 


columns 

1-2 

3-4 

5-12 

13-72 

data 

N1 

LIST 

YAM 

Hollerith  information 

formats 

12 

12 

E8.5 

10A6 

1.  N1  specifies  the  crystal  as  follows: 


Crystal  System  N1 


triclinic  1 

monoclinic 2 
orthorhombic  3 

tetragonal  4 

hexagonal  5 

cubic  6 


2.  B  LIST  =  0  or  blank,  the  complete  reflection  data  with  observed  and  calcu¬ 
lated  [(sin  6)/\'^  are  printed  out.  If  LIST  ^  0,  only  the  cell  parameter  and  errors 
are  printed. 

3.  YAM  is  the  wavelength  of  x  radiation  used,  or  else  YAM  is  0.  If  YAM  =  0 
or  blank,  the  program  assumes  that  (sin  e)/K  values  have  been  loaded.  If  YAM  =  \, 
the  program  assumes  that  sin  e  values  have  been  loaded. 

B.  Reflection  Card: 


columns 

1-4 

5-8 

9-12 

13-22 

23-32 

data 

h 

k 

1 

S 

W 

formats 

14 

14 

14 

E10.5 

E10.5 

1.  h,k,'tare  the  Miller  indices 

2.  S  is  either  sin  B  or  (sin  e)/\  (see  above) 

3.  w  is  a  least-squares  weighting  factor.  B  w  is  not  punched,  the  program 
assumes  w  »  (sin  b)/k. 

*lt  is  always  assumed  by  the  program  that  c  is  the  unique  axis;  thus  the  first  setting 
is  used  for  the  monoclinic  system. 
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The  reflection  card  deck  is  terminated  with  a  blank  card. 

A  second  calculation  may  follow.  The  reflection  deck  of  the  last  calculation  is 
terminated  with  two  blank  cards. 


1.3  Method  of  Calculation 

Let  S*  =  4(iin»fl)A* 

Hjj  =  hjhj,  i,  j  =  1,2,3 
gii  =  al*!,  i,  j  =  1,2,3 

where  h|  and  are  respectively  a  Miller  index  and  a  reciprocal  cell  edge  vector. 
Then 

where  the  repetition  of  a  subscript  as  a  superscript  in  the  same  term  of  an  eiqsression 
denotes  summation  over  that  index. 

The  least  squares  criterion  requires  that 

F  =  ”  “ini***”. 

whence 


or 


[**!.■  "u]  =  1-2.3. 


where  the  brackets  denote  summation  over  all  e:qperimental  observaticms,  and  1  is  a 
welj^ting  factor,  different  In  general  for  each  observation. 

After  collecting  terms,  this  becomes 

(  2  -  S„.)  (  2  -  « 1, )  [l  H„.  Hi  J  g*J  =  (  a  -  S„.)  [fS*  H„  J 

where  n,in,l,  j  s  1,2,3;  n  ^  n;  j  ^  i;  and  is  the  Kroneker  delta.  In  matrix  notation  this 
becomes 


where 


(N"«XU))  ■  (°("-))  "  (“<"■)) 

*’(n«Xii)  ~  ~  ®n»)  (2  “  ®lj)  [**!,■  ®ll] 

G  -  b”* 

(n«)  ^ 

V)  =  ‘ 

The  are  found  by  the  usual  matrix  methods.  The  variances  likewise  are  found  by  the 
usual  methods,  i.e., 

^0  =  -L  . 

("■)  |h|  k-< 
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where  is  the  cofactor  of  h/„a)(na)  element  of  |h|,  k  Is  the  number  of  terms  In  the 

bracket  summation,  and  <C  is  the  oraer  of  |h|. 

The  uncertainties  in  the  unit  cell  parameters  and  volume  are  then  found  by  error 
propagation  methods.  The  equations  are  quite  messy  and  will  not  be  reproduced  here. 

1.4  Output 

The  ou4>ut  is  sell  explanatory.  It  consists  of  the  reciprocal  cell  parameters,  the 
real  cell  parameters  with  uncertainties,  the  real  cell  volume  with  uncertainty,  and,  if 
desired,  the  observed  and  calculated  s  values  are  listed  for  each  reflection. 

1.5  Limitations 

No  more  than  ISO  observations  may  be  used. 

1.6  ^cial  Subroutines  Called 

In  addition  to  the  executive  program,  the  following  special  subroutines  are  used: 


a.  PARAM  c.  ERRREL 

b.  ERRPRM  d.  OUTPRM. 


Other  nonlibrary  programs  used  are: 

e.  RECIP 

f.  ARCSIN 

g.  Modified  versions  of  RW  MATS  and  RW  DET  (SHARE  distn.  no.  635). 


1.7  %)ecial  Tape  Requirements 
None 


1 


2.  QUASI  NORMAUZATION  OF  STRUCTURE  FACTORS 
2.1  Purpose 

Prepare  a  BCD  tape  containing  h,k,  I,  and  8^-1,  similar  to  the  one  written  by  XRDDR, 
using  for  input  the  h,k,'t,  tape  written  ^  XRDDR  in  the  Busing  least-squares  program 
format  (1). 


2.2  Input 

Ihe  input  consists  of  a  card  input,  which  is  loaded  following  the  program  deck,  and 
a  t«q>e  input,  which  usually  will  consist  of  the  data  tape  output  from  SUBROUTINE  OUTPUT 
of  XRDDR.  The  Uq>e  input  is  from  logical  tape  16. 
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A.  Control  Card: 


1.  NOIB  is  the  number  of  reflection  records  on  logical  tape  16. 

2.  NASP  is  the  number  of  different  atomic  species  In  the  crystal. 

3.  NF  is  the  number  of  files  on  the  input  tape  to  be  spaced  over.  This  allows 
input  to  be  taken  from  a  library  tape. 

4.  BTl  is  the  modified  Wilson  equation  temperature  factor  B  (as  taken  directly 
from  the  output  of  XRDDR). 

5.  EX  is  the  exponent  in  the  modified  Wilson  equation  (labeled  X  in  the  output  of 
XRDDR). 

6.  LIB  is  1  if  a  library  tape  is  used  as  input.  It  is  blank  otherwise,  (tt  is 
assumed  that  the  first  record  of  a  library  tape  is  an  identification  record.) 

Atomic  Scattering  Factor  Cards: 


columns 

1-8 

9-16 

17-24 

25-32 

33-40 

41-44 

data 

AS 

ASl 

BS 

BSl 

CS 

AN 

formats 

F8.4 

F8.4 

F6.4 

F8.4 

F8.4 

F4.0 

The  quantities  AS,  ASl,  BS,  BSl  and  CS  are  the  parameters  A,  a,  B,  b,  and  c 
defined  and  listed  by  Forsyth  and  Wells  lor  calculation  of  atomic  scattering  factors  (3). 
AN  is  the  number  of  atoms  of  the  atomic  species  in  the  unit  cell.  There  is  one  card 
for  every  different  atomic  species  in  the  crystal.  The  deck  consists  cf  one  control 
card  followed  by  NASP  atomic  scattering  factor  cards. 

B.  Reflection  Data  Tape: 


columns 

1-9 

10-18 

19-27 

28-36 

37-54 

55-63 

datA 

h 

k 

4 

F* 

blank 

(•in  6)/K 

formats 

F9.2 

F9.2 

F9.2 

F9.2 

18X 

F9.6 

2.3  Method  of  Calculation 

f2  values  are  converted  to  6^  values  according  to  the  following  equation 

Hie  reflections  are  read  from  the  input  tape  and  processed  in  groups  of  SOO  (or  less 
for  the  last  group).  Thus,  any  number  of  reflections  may  be  processed  in  one  computer 
run. 


2.4  Output 

Logical  tiqie  15  is*  rewound  at  the  beginning  of  the  calculation.  Output  is  BCD  on  t^;>e  15. 
The  program  writes  an  END  FILE  on  15  and  rewinds  it  when  the  calculation  is  finished. 


columns 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

data 

NOIB 

NASP 

NF 

BTl 

EX 

LIB 

formats 

110 

no 

no 

E10.4 

E10.4 

no 

II 
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columns 

1-4 

5-8 

9-12 

13-32 

data 

h 

k 

1 

£2-1 

formats 

14 

14 

14 

F20.5 

2.5  Limitations 

There  can  be  no  more  than  twenty  different  atomic  species  in  the  crystal. 

2.6  Social  Subroutines  Called 

A  FAP  coded  subroutine  FLSKPD  is  called  which  forward-spaces  tape  16  the  number 
of  files  specified  by  NF  (see  control  card). 

2.7  ^cial  Tape  Requirements 

Logical  tape  16  is  used  for  BCD  reflection  input. 

3.  RATIONAL  DEPENDENCE* 

3.1  Purpose 

To  determine  the  extent  of  rational  dependence,  as  defined  by  Hauptman  and  Karle  (4), 
in  a  crystal  structure  using  the  normalized  structure  factor  magnitudes. 

3.2  Input 

The  input  consists  of  a  card  input,  which  is  loaded  following  the  program  deck,  and 
a  tape  input,  which  usually  will  consist  of  the  data  tape  output  from  program  SF  NORM  of 
XRDDR.  The  tape  input  is  from  logical  tape  12. 

A.  Card  Input 

There  are  two  types  of  cards  in  the  card  input. 

a.  72  columns  of  Hollerith  characters  (these  characters  are  used  to  title  the 
output). 

b.  Control  Card; 


colunms 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 

data 

NA 

NINT 

NFIN 

MINN 

AN 

SDMIN 

JPUT 

formats 

no 

no 

no 

no 

E10.4 

E10.4 

no 

1.  NA  is  1  for  noncentrosymmetric  crystals;  it  is  0  or  blank  for  centrosymmetric 
crystals. 

2.  NINT  is  the  initial  modulas. 


♦After  Block  and  Yannoni,  National  Bureau  of  Standard  (unpublished). 
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3.  NFIN  is  the  final  modulas. 

4.  MINN  is  the  minimum  number  of  reflections  considered  significant  in  a  sub¬ 
set  showing  rational  dependence  effects.  Results  will  not  be  output  for  subsets  of 
reflections  with  fewer  than  MINN  members. 

5.  AN  is  the  minimum  number  of  standard  deviations  which  is  considered  a 
significant  difference  between  the  subset  average  and  overall  average  of 

6.  SDMIN  is  the  smallest  deviation  between  the  subset  average  and  the  overall 
average  of  E^,  where  the  deviation  is  considered  to  significantly  show  the  presence  of 
rational  dependence  effects. 

7.  JPUT  is  not  0  when  complete  listings  of  reflections  are  desired  for  each  sub¬ 
set  for  which  evidence  of  rational  dependence  has  been  found.  JPUT  is  0  when  the 
reflection  listings  are  not  desired. 

Additional  control  cards  may  follow  if  additional  calculations  are  desired.  The  card 
deck  is  terminated  with  a  blank  card. 

B.  Tape  Input. 

This  tape  may  be  tape  number  5  or  7  from  program  SF  NORM  of  XRODR: 

1.  The  first  record  on  the  tape  consists  of  the  number  of  reflection  records  on 
the  tape  (format  17). 

2.  Reflection  records: 


columns 

1-4 

5-8 

9-12 

10-32 

data 

h 

k 

E^-l 

formats 

14 

14 

14 

F20.5 

3.3  Method  of  Calculation 

A  FORTRAN  source  program  written  by  Block  and  Yannoni  was  obtained  from  Dr.  Block 
of  the  National  Bureau  of  Standards.  The  program  has  been  modified  with  respect  to  reflec¬ 
tion  input;  certain  redundancies  in  the  selection  of  subsets  have  been  eliminated;  a  pro¬ 
vision  for  specifying  the  minimum  significant  number  of  reflections  in  the  subsets  has 
been  added;  the  criterion  for  accepting  as  significant  a  deviation  of  average  E^  from  the 
overall  value  has  been  operationally  improved  by  making  it  independent  of  the  numbers  of 
reflections  in  the  subsets;  and,  finally,  the  output  has  been  vastly  e:qpanded  and  improved. 
Otherwise  the  calculations  are  done  essentially  as  coded  by  Block  and  Yannoni. 

All  reflections  in  a  subset  satisfy  the  relation 

ah  +  bk  +  C't  ■  n  (mod  m) 

where  m  is  a  positive  integer  limited  (by  the  program)  to  a  maximum  value  of  14;  a,  b, 
and  c  may  have  integral  values  ranging  from  -13  to  +14,  depending  upon  the  value  of  m. 

If  the  value  of  average  E^  for  a  subset  differs  significantly  from  the  average  over  all 
observed  E^,  then  rational  dependence  is  said  to  exist  in  the  crystal  structure. 

The  program  automatically  checks  all  possible  subsets  for  rational  dependence  effects 
for  all  moduli  (i.e.,  ^  m)  between  the  limits  specified  on  the  control  card. 
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3.4  Output 

The  output  consists  of  three  parts  A,  B,  and  C.  A  and  B  are  always  obtained,  and  C 
Is  or  Is  not  obtained  at  the  option  of  the  program  user. 

A.  This  is  a  leading  one-page  printout  consisting  of  the  control  information,  overall 
average  E^,  and  the  total  number  of  reflections  on  the  input  tape. 

B.  For  each  subset  for  which  significant  rational  dependence  effects  are  found  the 
following  is  printed:  a,b,c,n,ni,  average  E^,  and  the  number  of  reflections  in  the  subset. 

C.  B  the  control  datum  JPUT  is  not  zero,  data  for  all  reflections  in  the  subset  defined 
by  the  printout  of  B  above  are  listed,  immediately  following  the  B  printout,  in  the  form: 
h,k,^,  E^l,  and  E^-A,  where  A  is  the  overall  average  E^. 


3.S  Limitations 

The  only  limitation  in  addition  to  those  already  mentioned  in  section  3.3  is  that  the 
number  of  input  reflections  be  less  than  5000. 


3.6  Special  Subroutines  Called 

A  subroutine  OUTRAT,  especially  designed  to  write  output  for  this  program,  is  called. 

3.7  l^clal  Tape  Requirements 

Logical  tape  12  is  used  for  reflection  input. 


4.  SIGMA-2  LISTINGS 
4.1  Purpose 

For  each  member  of  an  ordered  list  of  reflections,  where  the  member  is  defined  by 
the  Miller  index  triple  h  and  norm^ized  structure  factor  magnitude  E,  the  program  lists 
all  pairs  of  reflections  h^,  £^,  and  h  j,  E^  from  a  similar  ordered  list  of  reflections  which 
satisfy  the  relation 

h  =  hj  +  h. . 

A  pair  of  reflections  and  hj  is  said  to  be  a  interaction  pair  for  the  reflection  h  if  the 
equation  above  is  satisfied. 


4.2  Input 

The  data  input  consists  of  one  control  card,  which  is  positioned  at  the  end  of  the  pro¬ 
gram  deck,  and  a  list  of  h  and  E^-1,  which  is  read  in  either  from  logical  tape  11  or  from 
cards.  The  ou^ut  t^s  from  program  SF  NORM  of  XRDDR  may  be  used. 


8 


NAVAL  KfSEARCH  LABORATORY 


A.  Control  Card: 


colunms 

1-10 

11-20 

21-30 

data 

Z 

IZ 

NOIB 

formats 

E10.4 

no 

no 

1.  All  reflections  for  which  E^>1  <  Z  are  rejected  from  consideration. 

2.  If  IZ  =  0  the  reflection  data  are  read  from  logical  tape  11.  If  IZ  0  the  reflec¬ 
tion  data  are  read  from  cards. 

3.  Lists  of  £,  Interaction  pairs  are  tabulated  for  the  NOIB  reflections  with  larg¬ 
est  E^-1  values.  Q  NOIB  Is  negative,  lists  are  prepared  for  the  entire  group  of  reflec¬ 
tions  for  which  E^-1  i  Z  (see  method  of  computation  below). 

B.  Reflection  Data. 

The  reflection  data  Is  prefaced  with  one  card  (or  record)  NRT  containing  the 
number  of  reflections  to  be  read.  The  format  Is  17. 

The  reflection  records  contain  the  following  information: 


columns 

1-4 

5-8 

9-12 

13-32 

data 

h 

k 

E»-l 

formats 

14 

14 

14 

F20.5 

4.3  Method  of  Calculation 

As  each  reflection  is  read  into  the  computer  it  Is  accepted  or  rejected  on  the  basis 
of  whether  or  not  the  criteria  on  h,k,'C  specified  In  REJVEC  (see  section  4.5)  are  satis¬ 
fied.  If  the  criteria  are  not  satisfied,  REJVEC  sets  E^-l  =  -10.  Then,  the  reflection  is 
tested  on  the  basis  of  the  magnitude  of  E^-l,  and  all  reflections  for  which  E^-l  <  Z  are 
rejected.  (The  minimum  physically  meaningful  value  of  E^.i  is  -l.) 

The  selected  set  of  reflections  Is  arranged  in  decreasing  order  of  magnitude  of  E^-1. 
Then  lists  of  Xj  interaction  pairs  are  calculated  and  tabulated  for  the  top  NOIB  reflections 
in  the  selected  set.  If  NOIB  has  been  entered  as  a  negative  number,  the  lists  are  calcu¬ 
lated  for  every  member  of  the  selected  set. 


4.4  Limitations 

The  selected  set  of  reflections  may  not  exceed  2400  reflections.  There  may  not  be 
more  than  1000  I,  interaction  pairs  for  each  h. 


4.5  l^cial  Subroutines  Called 
a.  REJVEC. 

This  subroutine  Imposes  acceptance  criteria  on  h,k,<(  of  each  reflection.  Hie 
program  user  may  write  his  ovm  p.'^gram  or  use  the  dummy  progrsun  already  written. 
The  quantities  Nl,  N2,  N3,  and  Q  in  the  calling  sequence  are  h,k,'(,  and  E^-1. 
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b.  SIGVEC. 

A  crystal-system-dependent  subroutine  which  selects  h.  and  hj  for  each  h.  If 
SIGVEC  for  the  trlclinic  case  Is  used,  and  the  reflection  decM  are  always  “blownup” 
so  as  to  include  all  symmetry  mates  of  all  reflections,  then  this  subroutine  becomes 
noncrystal-system  dependent.  SIGVEC  versions  for  triclinic,  monoclinic,  and  ortho¬ 
rhombic  crystals  are  available. 

c.  OUTSIG. 

l^is  subroutine  prints  out  the  list  of  h^,  andh^,  E.  and  the  products  E.E.Efor 
each  h,  E.  '  ’ 


4.6  fecial  Tape  Requirements 

Logical  tape  11  is  used  for  reflection  input. 

5.  TRIPLE  PRODUCT  SUMMATION  FOR  ORTHORHOMBIC  CRYSTALS 
5.1  Purpose 

These  programs  (one  for  centrosymmetric  and  one  for  noncentrosymmetric  crystals) 
are  designed  to  calculate  the  quantity  E.EjE,,  as  defined  by  Hauptman  and  Karle  (5),  or  the 
quantity  lEjEjEj  |  cos  (<^i  +  02+4i3)>  as  defined  by  Karle  and  Hauptman  (6). 


5.2  Input 

The  input  consists  of  a  card  input,  which  is  loaded  following  the.  program  deck,  and 
a  tape  input,  which  usually  will  consist  of  an  output  tape  from  program  SF  NORM  of 
XRDDR. 

A.  Control  Card; 


columns 

1-10 

11-20 

21-30 

31-40 

data 

NA 

IM 

JM 

KM 

formats 

no 

no 

no 

no 

1.  NA  is  the  number  of  atoms  (exclusive  of  hydrogen  atoms)  in  the  unit  cell. 

2.  IM  is  the  maximum  value  of  Miller  index  h. 

3.  JM  is  the  maximum  value  of  Miller  index  k. 

4.  KM  is  the  maximum  value  of  Miller  index  f . 

Triple  Product  Selection  Card: 


columns 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

data 

hi 

hi 

-tl 

^2 

*^2 

^2 

formats 

no 

no 

no 

no 

no 

no 
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1.  hj,ku>C|  are  the  Miller  indices  of  the  reflection  with  normalized  structure 
factor  E|  (see  section  S.l). 

2.  hjfkj,  <t2are  the  Miller  Indices  of  the  reflection  with  normalized  structure 
factor  E,  (see  section  5.1). 

The  deck  consists  of  one  control  card  followed  by  any  number  of  triple  product  selec¬ 
tion  cards.  Calculation  is  terminated  byjhe  second  blank  card  read  (the  first  blank 
card  causes  calculation  of  the  triple  for  hj.hj*  0). 

B.  Reflection  Data. 

The  first  record  on  the  reflection  data  input  tape  contains  the  number  of  reflec¬ 
tion  records  to  follow  (format  17). 

Each  reflection  record  consists  of: 


colunms 

1^4 

5-8 

9-12 

13-32 

data 

h 

k 

i 

E*-l 

formats 

14 

14 

14 

F20.5 

Tape  input  is  from  logical  tape  9. 

For  many  problems  it  may  be  necessary  to  recompile  the  source  program  in 
order  that  the  dimension  assignments  on  the  triply  indexed  variable  E  be  changed  (see 
section  5.3). 


5.3  Method  of  Calculation 

Hie  program  stores  each  value  of  (e^1)|,,^{  as  the  triply  indexed  FORTRAN  variable 
E(h,k,  1).  The  program  then  runs  through  the  list  of  reflections  in  a  triple  nest  of  DO 
loops  in  order  to  select  summands. 

The  orthorhombic  symmetry  operations  on  h,k,'(  are  coded  into  the  program. 

Since  the  program  uses  h,k,  i  as  Indices,  these  quantities  always  must  be  positive. 
The  program  automatically  rejects  one-  and  two-dimensional  data.  Core  size  limits  the 
maximum  values  of  h,k,  I  rather  severely.  As  listed  in  this  report,  the  program  accepts 
maximum  values  of  h,k,'{  as  large  as  30,  10,  30.  E  is  the  only  dimensioned  variable  in  the 
programs,  and  most  of  the  core  is  available  to  it;  however  it  will  frequently  be  necessary 
to  recompile. 

The  only  difference  between  the  programs  for  centrosymmetric  and  noncentrosym- 
metric  crystals  is  that  for  the  former,  the  average  is  divided  by  NA^V8«  whereas  the 
latter  is  divided  by  NA^V2.  (Note  the  third  statement  following  statement  number  8  in 
the  FORTRAN  listing  given  in  Appendix  A.) 

Whereas  only  three-dimensional  data  are  included  in  the  summations,  h| ,  ii],  >13  need 
not  be  three  dimensional. 

If  none  of  the  Miller  indices  of  hj,  ii.,  ii,  are  zero,  the  program  includes  the  correc¬ 
tion  term  (E*h,+  S*h,+  K*h,"  2)/nX‘^*  in  uie  results. 
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5.4  Output 

The  output  for  each  reflection  includes: 
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hj,  hi,hj 

b.  the  number  of  terms  in  the  sum, 

c.  the  scaled  triple  product  average  without  correction  term. 

d.  the  scaled  triple  product  average  with  correction  term. 

e.  the  correction  term. 

5.5  Limitations 
See  section  5.3. 

5.6  %)eclal  Subroutines  Called 
None. 

5.7  Social  Tape  Requirements 

Logical  tape  9  is  used  for  reflection  input. 

6.  STRUCTURE  FACTOR  CALCULATION 

6.1  Purpose 

To  provide  a  convenient  means  for  calculation  of  structure  factors  and  quasi- 
normalized  structure  factors  for  use  during  preliminary  stages  of  structure  determination. 
Reflection  input  is  from  the  E^-1  or  output  tapes  of  XRDDR.  Isotropic  temperature  fac¬ 
tors  are  used. 

6.2  Input 

The  input  may  be  divided  into  general  and  control  information,  and  reflection  data. 

A.  General  and  Control  Information. 

The  cards  in  this  part  are  included  with  the  program  deck. 

1.  72  columns  of  Hollerith  characters  (format  12A6).  These  Hollerith  characters 
are  used  to  title  the  output. 

2.  Control  Card  No.  1: 


columns 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 

IC 

IE 

NR 

NF 

IP 

LIB 

LINE 

formats 

no 

no 

no 

no 

no 

no 

no 
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IC  -  0  if  crystal  Is  centric;  IC  0  if  crystal  is  acentric. 

IE  s  0  if  F  3  reflection  data  input  is  used;  IE  ^  0  if  E^.!  reflection  data  input  is 
used. 

NR  is  the  number  of  input  reflections  (if  E^-l  reflection  data  is  used,  this  may 
be  left  blank). 

NF  is  the  number  of  files  to  be  spaced  over  if  the  reflection  data  are  read  from 
a  library  tape. 

IP  =  0:  E  output  mode  is  specified;  IP  /  0:  F  output  mode  is  specified.  IP  is 
irrelevant  for  centric  crystals  (see  section  6.5.B  for  a  description  of  the  modes). 

LIB  =  0;  library  tape  not  used  for  reflection  input;  LIB  0:  library  tape  is  used 
for  reflection  input. 

Line  =  0;  reflection  data  is  read  from  logical  tape  10;  LINE  ^  0:  reflection  data 
is  read  from  cards. 

3.  Cell  Constant  Card: 


columns 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

data 

a 

b 

c 

a 

y 

formats 

E10.4 

E10.4 

E10.4 

E10.4 

E10.4 

E10.4 

a,b,c,a,/3,r  Are  the  real-space  unit  cell  parameters  in  units  of  angstroms  and 
degrees. 

4.  Control  Card  No.  2: 


columns 

1-10 

11-20 

21-30 

31-40 

31-50 

data 

BT 

XT 

SFC 

KKT 

lOUT 

formats 

E10.4 

E10.4 

E10.4 

no 

no 

BT  and  XT  are  the  parameters  B  and  X  of  the  modified  Wilson  equation  as  taken 
from  the  output  of  XRDDR  (if  KKT  =  1,  these  fields  may  be  left  blank). 

SFC  is  a  scale  factor  by  which  all  Fg  arc  to  be  multiplied  to  bring  them  to  the 
proper  scale.  It  is  irrelevant  if  E^data  are  input. 

KKT  has  the  value  of  0,1,  or  2,  which  serves  to  select  the  desired  structure  fac¬ 
tor  normalization  procedure  as  follows: 

KKT  =  0:  apply  exact  normalization  procedure  to  F^;  apply  approximate  normali¬ 
zation  procedure  to  Fq  or  E^. 

KKT  =  1:  apply  exact  normalization  procedure  to  F^;  apply  exact  normalization 
procedure  to  Fg  or  Eg. 

KKT  =  2:  apply  approximate  normalization  procedure  to  F^;  apply  approximate 
normalization  procedure  to  Fg  or  Eg. 

(see  section  6.4  for  explanations  of  exact  and  approximate  normalization.) 

If  lOUT  ^  0,  a  BCD  output  (logical  tape  11)  will  be  written  (see  section  6.6). 
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5.  Scale  factor  cards. 

H  IE  is  0  (i.e.,  Fg  reflection  data  are  input),  a  deck  of  scale  factor  cards  is  input. 
Each  card  contains  one  scale  factor  (format  E10.4).  The  deck  is  terminated  with  a 
blank  card.  The  scale  factors  are  Indexed  in  the  order  loaded  and  are  applied  accord¬ 
ing  to  the  procedure  given  in  section  6.2.B  under  description  of  SCF. 

6.  Atomic  scattering  factor  cards. 

There  is  one  card  for  each  unique  atomic  species  in  the  crystal: 


columns 

1-8 

9-16 

17-24 

25-32 

33-40 

data 

A 

a 

B 

b 

C 

formats 

F8.4 

F8.4 

F8.4 

F8.4 

F8.4 

A,a,B,b,  and  c  are  the  parameters  used  in  the  Vand,  Eiland,  Pepinsky  approxi¬ 
mation  equation  for  calculation  of  atomic  scattering  factors  as  defined  by  Forsyth 
and  Wells  (3,7). 

The  deck  of  atomic  scattering  factor  cards  is  terminated  by  a  blank  card. 

The  loading  order  of  the  atomic  scattering  factor  cards  is  important  and  must 
be  consistent  with  the  quantities  I  on  the  coordinate  cards  (see  section  6.2.A.7). 

7.  Coordinate  cards. 

There  is  one  card  for  each  atom  in  the  asymmetric  unit: 


columns 

1-10 

11-20 

21-30 

31-40 

41-50 

data 

X 

Y 

Z 

B 

I 

formats 

E10.4 

E10.4 

E10.4 

E10.4 

no 

X,Y,  and  Z  are  the  fractional  atomic  coordinates  of  one  atom."' 

B  Is  the  isotropic  temperature  factor. 

I  is  an  integer  which  relates  the  atom  to  one  of  the  atomic  scattering  factor 
cards  and  thus  specifies  the  atomic  type.  Each  atomic  scattering  factor  card  is 
numbered  1,2,3,...  in  sequence  as  it  is  read  into  the  computer.  Thus,  if  the  second 
scattering  factor  card  to  be  read  is  for  carbon,  then  I  =  2  on  all  carbon  coordinate 
cards. 

The  deck  of  coordinate  cards  is  terminated  with  a  blank  card. 

8.  Transformation  cards. 

The  atoms  in  the  asymmetric  unit  are  transformed  according  to  the  relations 
given  in  *International  Tables  for  X-Ray  Crystallography,”  Vol.I,  (N.F.M.  Henry  and 
K.  Lonsdale,  Birming^m^Kynoch  Press,  1952)  so  as  to  fill  out  the  unit  cell.  There 
is  one  card  for  each  transformation. 


columns 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

data 

T1 

T2 

T3 

U1 

U2 

U3 

formats 

E10.4 

E10.4 

E10.4 

E10.4 

E10.4 

E10.4 

'*'If  an  atom  is  at  the  origin,  X,Y,Z  should  be  given  the  values  1,1,1  instead  of  0,0,0. 
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The  transformations  are  done  as  follows: 


X' 

»  T1 

+  Ul-X 

Y' 

II 

+  U2-Y 

V 

=  T3 

+  U3-Z. 

Example:  %)ace  group  P2|/c 

Transformations:  x,y,(;  x,l/2 

■  y» 

1/2  +  I ;  X, 

1/2  +  y,  1/2 

Transformation  cards 

T1  T2 

T3 

U1  U2 

U3 

.0  .0 

.0 

-1.  -1. 

-1. 

.0  .5 

.5 

1.  -1. 

1. 

.0  .5 

.5 

-1.  1. 

-1. 

The  deck  of  transformation  cards  is  terminated  with  a  blank  card. 

B.  Reflection  Data. 

The  reflection  input  may  be  Included  along  with  the  program  deck,  following  the 
general  and  control  Information  cards,  or  be  taken  from  logical  tape  10.  It  may  be 
taken  directly  from  an  XRDDR  output  tape  or  from  a  library  tape,  each  file  of  which 
contains  the  data  from  an  XRDDR  input  or  output  tape.  In  the  latter  case,  the  pro* 
gram  assumes  that  each  fUe  is  prefixed  with  one  identification  record. 

The  program  calculates  structure  factors  for  each  reflection  in  the  reflection 
input. 

If  reflection  input  is  used  (tin  9)/k  values  may  or  may  not  be  included  (they 
are  automatically  Included  on  the  X^DR  output  Uqpes).  The  format  is  the  same 
as  that  used  lor  input  to  the  Busing  least -squares  program. 

S  E^-1  reflection  input  is  used,  the  first  record  must  consist  of  the  number  of 
reflections  to  be  read  (format  17).  No  provision  is  made  for  input  of  (tin  0)/k  data. 

reflection  data: 


columns 

1-9 

10-18 

19-27 

28-36 

37-45 

46-54 

55-64 

h 

k 

1 

Fo 

(blank) 

SFC 

(tin  S)/X 

formats 

F9.2 

F9.2 

F9.2 

F9.2 

9X 

F9.2 

F9.6 

The  quantity  SCF  is  a  floating  point  integer  which  relates  the  reflection  to  one  of 
the  scale  factors  described  in  section  6.2.A.5  above.  The  integer  corresponds  to  the 
loading  sequence  number  of  the  scale  factor  (i.e.,  the  first  scale  factor  loaded  corre¬ 
sponds  to  SCF  s  1.0,  etc.).  SCF  may  be  punched  on  every  card.  Alternatively,  the 
reflection  input  may  be  grouped  into  sets  where  all  reflections  in  a  set  are  multiplied 
by  the  same  scale  factor.  In  this  case,  it  is  sufficient  to  punch  SCF  in  the  first  card 
only  of  each  set.  This  is  the  same  scsiing  procedure  used  in  the  Busing  least-sqLiares 
refinement  program;  however,  in  this  program,  the  Fg  values  are  scaled. 
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e‘-1  reflection  data: 


columns 

1-4 

5*8 

9-12 

13-32 

Hata 

h 

k 

Ko*-l 

formats 

14 

14 

14 

F20.5 

Jt  another  calculation  Is  to  follow,  a  card  with  a  1  punched  In  colunm  7  follows 
the  last  card  of  the  main  deck;  otherwise,  a  blank  cstrd  Is  last. 

6.3  Data  Deck 

The  composition  of  the  data  deck  Is  as  follows: 

a.  Hollerith  card. 

b.  control  card  no.  1. 

c.  cell  constant  card. 

d.  control  card  no.  2. 

e.  all  scale  factor  cards  (if  IE  ^  0,  there  are  no  scale  factor  cards). 

f.  blank  card. 

g.  all  atomic  scattering  factor  cards. 

h.  blank  card. 

1.  all  coordinate  cards. 

j.  blank  card. 

k.  all  transformation  cards  (If  any). 

l.  blank  card. 

m.  reflection  cards  unless  reflection  input  is  from  tape. 

n.  a  card  with  1  in  column  7  or  a  blank  card  depending  upon  whether  or  not 
another  calculation  Is  to  follow. 


6.4  Method  of  Calculation 

Structure  factors  are  calculated  from  the  total  unit  cell  contents  (rather  than  from 
the  asymmetric  unit  contents)  using  the  equations  for  the  triclinic  case.  Nevertheless, 
coordinates  may  be  entered  for  one  asynunetrlc  unit  only,  provided  that  transformations 
follow  by  which  the  asymmetric  unit  set  is  e;qMUKled  to  fill  the  unit  cell.  The  transforma¬ 
tions  may  be  celled  directly  from  the  *Intemational  Tables  for  Crystallography,*  Vol.  I 
(1952).  hi  this  way,  the  program  is  completely  general  in  its  applicability  while  requiring 
the  most  sinqile  t^  of  input. 

The  quasi  imrmalization  of  structure  factors  (and  in  the  case  where  E^-l  reflection 
input  is  used,  the  inverse  normalization)  may  be  done  in  either  of  two  ways,  designated 
as  exact  and  approximate  normalization. 


1.  Exact  normalization: 
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2.  Approximate  normalization: 

j  F*  exp  [BT(iin2 

where 

fi  =  fo.i  exp[-Bj(sin2  ^)A*] 
and  N  is  the  number  of  atoms  in  the  unit  cell. 


6.5  Printed  Output 

A.  Centric  Crystals. 

The  output  for  each  reflection  is: 

a.  h,k,'C 

b.  Fo 

d.  E, 

e.  E. 

f.  AF 

g.  AE. 

B.  Acentric  Crystals. 

Either  one  of  two  modes  of  output  may  be  selected: 
1.  F  mode 


a.  hiky-C 

b.  Eo.  Ec 

^0*  ®0 

e.  A^,B^ 

f.  hA,  AB 

g.  AF. 

2.  E  mode 


a. 

h,k,-t 

b. 

c. 

d. 

Co»  Cel 

e. 

Po.  be  J 

£. 

AC,  AD 

g- 

AE. 

where  1e|*=  C*+  d* 
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In  addition  to  the  reflection  output,  the  program  prints  £|AFl,  iIFqI,  and 
R  s  i;  I AF I  / 1 1  Fg  I  both  including  and  excludii^  unobserved  data. 


6.6  Tape  Output 

K  lOUT  4  0  the  program  writes  a  reflection  output  tape  containing  one  of  the  follow¬ 
ing  sets  of  data  for  each  reflection  (format  7F10.4): 

a.  Centrosymmetric  crystals 

h.k.^.Fg.Eg.AF.AE 

b.  Noncentrosymmetric  crystals 

1.  1P  =  0: 

h.k.-t.Cg.Dg.AC.AD 

2.  IP?^0: 

h.k,-l,Ag,Bg.AA,AB 

where  the  observed  quantities  are  scaled  and  given  the  signs  of  the  corresponding  calcu¬ 
lated  quantities. 


6.7  Social  Subroutines  Called 

a.  OUTSFC. 

An  output  routine  designed  for  these  calculations  only. 

b.  RECIP. 

c.  The  FAP  coded  tape  utility  program  written  by  P.  Gum.* 

d.  INTO. 

6.8  l^cial  Tape  Requirements 

a.  Logical  10  is  used  for  reflection  input. 

b.  Logical  11  is  used  for  reflection  output. 

7.  INTERATOMIC  DISTANCES  AND  ANGLES  WITH 
INTERPOLATION  FOR  PEAK  CENTER  LOCATION 

7.1  Puipose 

Given  a  basic  set  of  atomic  or  peak  coordinates  and  a  set  of  symmetry  transformations 
on  these  coordinates,  calculate  and  list  all  Interatomic  distances  b,  where  0.7  <  b  <  BMX, 
and  calculate  all  angles  between  connected  pairs  of  distances  for  all  possible  combinations 
in  any  crystal. 


■*‘Diffraction  Branch,  Optics  Division,  NRL. 
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7.2  Input 

The  data  input  is  from  cards  which  are  loaded  immediately  following  the  program 
deck.  The  cards  are  described  in  loading  sequence  as  follows: 

A.  72  columns  of  Hollerith  characters  (format  12A6).  These  Hollerith  characters  are 
used  to  title  the  output. 

B.  Cell  Constant  Card: 


columns 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

data 

a 

b 

c 

a 

y 

formats 

E10.4 

E10.4 

E10.4 

E10.4 

B10.4 

E10.4 

The  quantities  are  the  real-space  unit  cell  parameters  expressed  in 

units  of  angstroms  and  degrees. 

C.  Bond  Length  Limit  Card: 


columns 

1-10 

11-20 

data 

BMX 

AMX 

formats 

E10.4 

E10.4 

BMX. 

Interatomic  distances  b  for  which  0.7  ^  b  $  BMX  are  considered  acceptable  and 
are  tabulated  for  printing. 

AMX. 

Angles  for  all  pairs  of  acceptable  distances  sharing  a  common  atom,  for  which 
the  distances  of  both  members  of  the  pair  are  less  than  or  equal  to  AMX,  are  tabu¬ 
lated  for  printing. 

[WAfUHNG:  When  a  large  number  of  atoms  or  peaks  are  involved  and  BMX  and 
AMX  are  both  large  (i.e.,  greater  than  about  3.0A),  the  number  of  angles  printed  can 
be  very  large.] 

D.  Grid  Interval  Card: 

H  peak  maxima  coordinates  are  to  be  found  by  interpolation  between  grid  points 
of  a  Fourier  calculation  result,  this  card  contains  the  fractional  grid  increments  in 
the  directions  of  X,  Y,  and  Z. 


columns 

1-10 

11-20 

21-30 

xmc 

Ymc 

ZINC 

formats 

E10.4 

E10.4 

E10.4 

For  example,  if  the  Fourier  calculation  was  done  in  increments  of  1/30,  1/60, 
and  1/120  along  X,  Y,  and  Z,  respectively,  then 

Xmc  »  .03333,  YINC  -  .01666,  and  ZINC  >  .008333. 

If  interpolation  is  not  to  be  done,  this  card  must  be  blank. 


1 
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E.  Coordinate  Card: 


columns 

1-10 

11-20 

21-30 

data 

X 

Y 

Z 

formats 

E10.4 

E10.4 

E10.4 

X,  Y,  and  Z  are  the  fractional  coordinates  of  an  atom  or  peak.  There  is  one 
card  for  each  atom  or  peak  in  the  basic  (untransformed)  set.  B  interpolation  is  to  be 
done,  each  coordinate  card  is  followed  by  an  interpolation  card  (see  section  7.2.  F 
below).  If  interpolation  is  not  to  be  done,  the  interpolation  cards  are  absent  from 
the  deck. 

The  coordinate  card  deck  is  terminated  by  a  blank  card. 

F.  Interpolation  Card. 

Interpolation  cards  are  included  in  the  coordinate  deck,  one  following  each  coor¬ 
dinate  card,  if  and  only  if  card  D  is  not  blank  (see  section  7.2.D).  .2.4). 

The  fractional  coordinates  of  the  grid  point  closest  to  a  peak  maximum  in  a 
Fourier  calculation  result  are  given  on  the  coordinate  card.  Relative  to  th<M  central 
grid  point,  designate  the  Fourier  summation  value  at  nearby  grid  points  by  P(i,  j,k), 
where  i,j,  and  k  have  integral  values  representing  the  (signed)  n*  -abers  of  grid  inter¬ 
val  separations  of  the  point  from  the  central  one  in  the  directions  of  X,  Y,  and  Z, 
respectively. 


The  interpolation  card  contains  the  following: 


columns 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 

P(0,0,0) 

P(*1,0,0) 

P(1,0,0) 

P(0,-1,0) 

P(0,1,0) 

P(0,0,-1) 

P(0,0,1) 

formats 

B10.4 

E10.4 

E10.4 

E10.4 

E10.4 

E10.4 

E10.4 

If  the  user  wishes  to  interpolate  some,  but  not  all,  points,  then  values  of  P(0,0,0)  = 
100.0  and  all  other  P  =  1.0  may  bo  entered  for  those  points  that  are  not  to  be 
inteipolated. 

Note  that  P(0,0,0)  must  be  the  largest  P  on  the  card.  Also,  all  P  must  be  greater 
than  zero.  Negative  or  zero  P  can  be  accommodated  by  adding  a  constant  Increment 
to  all  P  such  that  the  results  are  all  positive. 

G.  Transformation  Card: 


columns 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

data 

T1 

T2 

TS 

U1 

U2 

U3 

formats 

E10.4 

E10.4 

E10.4 

B10.4 

E10.4 

E10.4 

There  is  one  transformation  card  for  every  symmetry  transformation  to  bo  applied 
to  the  basic  set  of  coordinates.  The  transformations  are  applied  according  to  the 
equations: 

X‘  e  T1  +  Ul-X 
Y‘  -  T2  +  U2*Y 
Z‘  -  T3  +  U3*Z. 

The  transformation  card  is  terminated  by  a  Wnnk  card. 
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7.3  Data  Deck 

The  composition  of  the  data  deck  is  as  follows: 

a.  Hollerith  card. 

b.  cell  constant  card. 

c.  bond  length  limit  card. 

d.  grid  interval  card,  or  blank  if  interpolation  is  not  desired. 

e.  all  coordinate  and  interpolation  cards. 

f.  blank  card. 

g.  all  transformation  cards  (if  any). 

h.  blank  card. 


7.4  Method  of  Calculation 


The  basic  set  of  coordinates  are  loaded,  Interpolated  (if  desired),  and  then  each  trans¬ 
formation  is  applied  to  the  basic  set,  generating  a  transformed  set  for  each  transforma¬ 
tion.  Thus,  if  there  are  N  atoms  in  the  basic  set  and  there  are  M  transformations,  the 
complete  set  contains  N(M  +  1)  atoms. 

Interatomic  distances  less  than  or  equal  to  BMX  are  tabulated  for  all  pairs  of  atoms 
in  the  basic  set  and  for  all  distances  ( <  BMX)  which  cross  the  transformed  set  boundaries 
for  the  complete  set. 

Atom  pairs  closer  together  than  0.7A  are  listed  in  a  separate  output  and  are  rejected 
from  further  consideration. 

Next,  the  list  of  accepted  interatomic  distances  is  scanned  to  find  pairs  of  interatomic 
distances  which  share  a  common  atom.  Angles  for  all  such  pairs  of  distances  are  tabu¬ 
lated  if  both  distances  are  less  than  or  equal  to  AMX.  Logical  tape  9  is  used  for  inter¬ 
mediate  storage  of  angles. 


Interpolation  is  done  by  fitting  the  seven  Fourier  grid  points  closest  to  the  peak  maxi¬ 
mum  to  a  Gaussian  function  and  then  finding  the  maximum  of  the  function.  The  corrections, 
£iX,  AY,  and  aZ  to  be  added  to  X,  Y,  and  Z  are  found  by  solving  the  matrix 


where 


etc. 


a'a  a*b  a-c 
a-b  b-b  b'c 


La-c  b’C  C'CJ 


lil 

P( 


LV0:.0)l/riog 


P(0,0,0) 


+  log 


P(0.0,0)‘ 

P(-1.0,0)J’ 


7.5  CXitput 

In  the  output,  each  atom  is  represented  by  two  Integers  separated  by  a  hyidien.  The 
first  Integer  specifies.the  transformation  (numbered  according  to  loading  order)  which 
has  been  applied,  and  the  second  Integer  specifies  the  atom  in  the  basic  set  (numbered 
according  to  loading  order)  from  which  tt  was  derived. 
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The  output  Is  in  the  following  order; 

a.  Atomic  coordinates  lor  the  complete  set. 

b.  Interatomic  distances. 

c.  Angles. 

d.  Interatomic  distances  less  than  0.7A. 


7.6  Limitations 

a.  The  complete  set  of  atoms  or  peaks  is  limited  to  500  in  number. 

b.  If  more  than  2900  acceptable  interatomic  distances  are  found,  BMX  is  decremented 
by  0.25A  and  the  calculation  is  begun  again.  This  process  is  repeated  until  a  total  of  2900, 
or  less,  acceptable  distances  are  found,  or  until  BMX  $  2.0.  In  the  latter  case,  distances 
and  angles  are  calculated  for  the  basic  set  of  atoms  only.  If  more  than  2900  acceptable 
distances  are  foimd  in  the  basic  set,  the  calculation  is  terminated  without  output  of  results. 

c.  No  more  than  200  distances  less  than  0.7A  can  be  accommodated. 

d.  There  is  no  limit  on  the  number  of  bond  angles. 


7.7  fecial  Subroutines  Called 

a.  subroutine  INTERP.  This  subroutine  does  the  interpolation  by  Gaussian  curve 
fitting. 

b.  subroutine  MATS  (SHARE  distr.  no.  635). 

c.  soubroutine  OUTBND:  (distance  and  angle  output). 

d.  function  DOTPRD. 

e.  function  ARCSIN. 

7.8  Social  Tape  Requirements 

Logical  Tape  9  is  used  for  intermediate  storage. 

8.  LEAST-SQUARES  PLANE  AND  LINE  FITTER 
8.1  Purpose 

This  program  fits  planes  and/or  lines  to  sets  of  points  and  calculates  angles  between 
planes  and  lines  which  have  been  fitted. 


8.2  biput 

a.  72  columns  of  Hollerith  characters  (format  12A6).  These  Hollerith  characters 
are  used  to  title  the  ou^t. 
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b.  Unit  Cell  Card: 


columns  1-10  11-20  21-30  31-40 


daU  A  b  c  a  /3 


formats  E10.4  E10.4  E10.4  E10.4  E10.4  I  E10.4 


The  quantities  aibiCiaj/Siy  are  the  real-space  unit  cell  parameters  e^qpressed 
in  units  of  angstroms  and  degrees. 

c.  Plane  or  Line  Reification  Card: 


columns 

10-3 

4-6 

1  1 

13-72 

data 

N 

NO 

mm 

HoUerith 

formats 

13 

13 

16 

10A6 

N  is  the  number  of  points  to  be  fitted. 

NO  is  the  plane  or  line  number. 

L  =  0  or  blank  if  the  points  are  to  be  fitted  to  a  plane. 

L  =  1  if  the  points  are  to  be  fitted  to  a  line. 

The  Hollerith  characters  in  colunms  13-72  are  used  as  a  page  heading  for  the 
line  or  plane  in  the  output. 

d.  Coordinate  Card: 


I-IO 

1 

El 

-30 

X 

z 

W 

E10.4 

1  E10.4  1 

E10.4 

E10.4 

X,  Y,  and  Z  are  fractional  coordinates  of  a  point  to  be  fitted  to  a  line  or  plane. 

W  is  a  least-squares  weighing  factor.  H  W  is  not  punched,  it  is  assumed  to 
have  a  value  of  one  by  the  program. 

e.  Angle  specification  card. 


This  card  specified  a  plane-plane,  line-line,  plane-line,  or  line-plane  pair  for 
which  the  dihedral  angle  is  to  be  calculated. 


NOl  is  the  number  of  the  first  plane  or  line. 

LI  s  0  or  blank  if  NOl  refers  to  a  plane. 

LI  -  1  if  NOl  refers  to  a  line. 

N02  is  the  number  of  the  second  plane  or  line. 
L2  -  0  or  blank  if  N02  refers  to  a  plane. 

L2  B  1  if  N02  refers  to  a  line. 
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The  plane  and  line  numbers  NO  are  used  to  identify  the  different  planes  and  lines 
both  internally  and  in  the  output.  The  numbering  may  be  different  from  the  order  of 
input.  No  two  planes  may  have  the  same  number,  and  no  two  lines  may  have  the  same 
number,  but  the  numbering  of  lines  and  planes  is  independent.  The  number  of  any 
plane  or  line  must  not  be  greater  than  SO. 


8.3  Data  Deck 

Ihe  composition  of  the  data  deck  is  as  follows: 

a.  Hollerith  card. 

b.  unit  cell  card. 

c.  plane  or  line  specification  card. 

d.  all  coordinate  cards  for  this  plane  or  line  (the  number  of  coordinate  cards 
must  be  the  same  as  the  N  punched  on  the  preceding  plane  or  line  specifi¬ 
cation  card). 

e.  plane  or  line  specification  card. 

f.  all  coordinate  caids  for  this  plane  or  line. 

g.  additional  sets  of  specification  and  coordinate  cards  for  as  many  planes  and 
lines  as  are  to  be  fitted. 

h.  blank  card. 

1.  all  angle  specification  cards,  if  any. 

j.  blank  card. 


8.4  Method  of  Calculation 

Ihe  method  described  in  detail  by  Shomaker,  Waser,  Marsh,  and  Bergman  (8)  is  used. 

The  program  is  written  in  the  fullest  generality  so  as  to  be  able  to  handle  points 
defined  in  any  three-dimensional  coordinate  system.  However,  no  provision  is  made  for 
symmetry  transformation  of  point  positi(His,  so  all  points  must  be  properly  transformed 
before  ii^t. 

For  a  given  coordinate  system,  the  program  sequentially  fits  all  desired  planes  and 
lines.  It  stores  the  appropriate  unit  vectors  for  the  planes  and  lines  by  plane  and  line 
number.  The  program  can  then  calculate  dihedral  angles  between  all  pairs  of  planes, 
lines,  or  plane-line  combinations  that  are  specified. 


8.5  Output 

The  ou^ut  is  largely  self  e]q>lanatory.  %)ecifically,  the  following  information  is 
given: 

a.  equations  of  the  planes  and  lines. 

b.  for  planes,  the  coordinates  of  the  centroid. 

c.  perpendicular  distances  of  each  point  from  the  least-squares  planes  and  lines. 

d.  average  deviations,  standard  deviations,  etc. 
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e.  cosines  of  angles  between  plane-plane  and  line-line  pairs. 

f.  sines  of  angles  between  line-plane  pairs. 

Each  point  is  numbered  in  the  output  according  to  loading  order.  Thus,  the  Hollerith 
section  of  the  plane  or  line  specification  card  (see  section  8.2.c)  may  be  used,  in  part,  to 
relate  the  points  to  atom  numbers  or  other  external  designations. 

8.6  Limitations 

a.  No  more  than  50  points  may  be  fitted  to  any  plane  or  line. 

b.  No  less  than  three  points  may  be  fitted  to  a  plane. 

c.  No  less  than  two  points  may  be  fitted  to  a  line. 

d.  There  may  be  no  more  than  50  planes  and  no  more  than  50  lines. 

e.  The  number  of  a  line  or  plane  may  not  be  larger  than  50. 

8.7  Special  Subroutines  Called 

a.  subroutine  RECIP. 

b.  subroutine  MTXMUL. 

c.  function  DOTPRD. 

8.8  %)ecial  Tape  Requirements 
None. 

9.  POINT  TO  PEAK  DISTANCE  CALCULATION 

9.1  Purpose 

This  program  provides  a  means  for  calculating  distances  from  one  or  several  points 
to  a  set  of  points  in  any  three-dimensional  coordinate  system. 

9.2  Input 

a.  72  columns  of  Hollerith  characters  (format  12A6).  These  characters  are  used  to 
title  the  output. 

b.  Unit  CeU  Card: 


The  quantities  are  the  real-space  unit  cell  parameters  in  units  of 

angstroms  and  degrees. 
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c.  Maximum  distance  card. 

This  card  contains  the  single  quantity  BMX.  All  distances  d  are  tabulated  for 
which  d  .$  BMX  (format  E10.4). 

d.  Point  Card: 


columns 

1-10 

11-20 

21-30 

data 

XC 

YC 

ZC 

formats 

E10.4 

E10.4 

E10.4 

XC,  YC,  and  ZC  are  the  fractional  coordinates  of  the  point  for  which  distances 
to  the  set  of  points,  specified  by  the  coordinate  cards  and  their  transformations,  are 
to  be  calculated. 

e.  Coordinate  card 

These  cards  contain  the  fractional  coordinates  X,  Y,  and  Z  of  the  basic  set  of 
points,  which  gives  rise  to  a  complete  set  via  the  transformations  (see  card  below), 
whose  distances  to  point  XC,  YC,  ZC  are  desired.  The  card  layout  is  the  same  as 
that  already  given  for  the  point  card  above. 

f.  Transformation  Cards: 


colunms 

1-10 

11-30 

21-30 

31-40 

41-50 

51-60 

data 

T1 

T2 

T3 

U1 

U2 

U3 

formats 

E10.4 

E10.4 

E10.4 

E10.4 

E10.4 

E10.4 

The  function  of  these  quantities  is  explained  in  section  7.2.G. 

9.3  Data  Deck 

a.  Hollerith  card. 

b.  unit  cell  card. 

c.  maximum  distance  card. 

d.  point  card. 

e.  all  coordinate  cards. 

f.  blank  card. 

g.  all  transformation  cards. 

h.  blank  card. 

i.  additional  point  cards  for  as  many  calculations  as  desired. 

j.  blank  card. 

9.4  Method  of  Calculation 

This  program  is  an  abbreviation  of  the  Interatomic  Distance  and  Angle  Program  (sec¬ 
tion  7).  It  does  not  provide  for  interpolation  or  angle  calculation. 
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f  9.5  Output 

The  output  of  distances  is  essentially  the  same  as  that  described  for  the  Interatomic 
Distance  and  Angle  Program  (section  7.5). 

9.6  Limitations 

a.  The  complete  set  of  points  must  not  exceed  500  in  number. 

b.  A  maximum  of  5000  distances  can  be  accommodated. 

9.7  fecial  Subroutines  Called 

a.  DOTPRD. 

9.8  Social  Tape  Requirements 
None. 

10.  FORM  FACTORS  FOR  THE  BUSING  LEAST-SQUARES 
REFINEMENT  PROGRAM 

10.1  Purpose 

Provide  an  automatic  process  for  preparing  form  factor  cards  used  for  input  to  the 
Busing  least-squares  refinement  program  ORXLS  (1). 

10.2  Input 

a.  One  blank  card. 

b.  Form  Factor  Parameter  Card; 


columns 

1-8 

9-16 

17-24 

25-32 

41-46 

data 

D 

a 

B 

b 

c 

HoUerith 

formats 

F8.5 

F8.5 

F8.5 

F8.5 

F8.5 

A6 

The  quantities  A,  a,  B,  b,  C  are  the  parameters,  as  defined  by  Forsyth  and 
Wells  (3),  used  in  calculating  the  form  factors.  The  six  columns  of  Hollerith 
characters  are  used  as  identification  on  the  output  cards. 

10.3  Data  Deck 

The  composition  of  the  data  deck  is  as  follows: 

a.  blank  card. 

b.  as  many  form  factor  cards  as  desired. 

c.  blank  card. 
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10.4  Method  of  Calculation 

The  analytical  approximation  of  Vand,  EUand,  and  Pepinsky  (7),  Is  used  as  extended 
by  the  work  of  Forsyth  and  Wells  (3). 

The  approximations  are  good  only  as  far  as  (sin  B)/k  1.40  (the  range  of  M^Ka  radiation). 
These  approximations  are  not*  recommended  for  highly  accurate  work. 

10.5  Output 

The  form  factor  card  images  are  written  on  logical  tape  10  for  peripheral  punching. 

For  each  atomic  species,  five  cards  are  punched  as  follows: 

a.  a  blank  separator  card. 

b.  four  cards  containing  the  form  factor  values  in  intervals  of  0.05  in  (sin  e)/\ 
from  1.55  to  0,  in  that  order. 

Each  card  has  an  F  in  colunm  72.  Columns  73-78  contain  the  Hollerith  characters 
specified  in  the  form  factor  parameter  cards.  Colunm  80  contains  the  sequence  number 
of  the  card,  l.e.,  1,2,3,  or  4. 

Essentially  the  same  information  also  is  printed. 


10.6  Special  Subroutines  Called 
None. 


10.7  ^cial  Tape  Requirements 

a.  Logical  tape  10  for  peripheral  punching. 

11.  VARIANCE-COVARIANCE  MATRIX  AND  ATOMIC  COORDINATE 
INPUT  FOR  THE  BUSING  FUNCTION  AND  ERROR  PROGRAM 

11.1  Purpose 

Ihis  program  produces  two  BCD  card  decks  containing,  respectively: 

a.  The  variance -covariance  matrix  for  a  set  of  atomic  position  coordinates, 
with  all  covariance  elements  given  a  value  of  zero. 

b.  Atomic  coordinates. 

These  cards  are  ready  for  input  to  the  Busing  Function  and  Error  Program  ORXFE  (2). 


11.2  Input 

a.  Preliminary  Card: 


columns 

1-9 

10-19 

data 

N 

SF 

formats 

19 

E10.5 

\ 


i 


! 


j 

'i 

\ 

I 

! 
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N  is  the  order  of  the  variance -covariance  matrix,  i.e.,  N  =  3x  (number  of  atoms). 

SF  is  a  scale  factor  by  which  the  matrix  is  to  be  multiplied.  Specifically, 

SF  =  I  *  (Fp  -  F^)V(m  -  n),  or  unity. 

b.  Unit  Cell  Card: 


columns 

1-10 

11-20 

21-30 

data 

a 

b 

c 

formats 

ElO.S 

ElO.S 

ElO.S 

If  the  standard  deviations  on  the  cards  described  below  are  in  fractional 
(dimensionless)  form,  a,b,  and  c  all  have  values  of  one. 

If  the  standard  deviations  are  in  units  of  angstroms,  then  a,b,  and  c  are 
the  unit  cell  edge  lengths  in  units  of  angstroms. 

c.  Standard  Deviation  Card: 


columns 

1-10 

11-20 

21-30 

data 

a(X) 

<^(Y) 

<^(Z) 

formats 

ElO.S 

ElO.S 

ElO.S 

The  quantities  <r(X),  (r(Y),  and  a(Z)  are  standard  deviations  or  uncertainties  of 
the  atomic  coordinates  X,  Y,  and  Z.  These  quantities  may  be  dimensionless  or 
may  have  units  of  angstroms  (see  section  11.2.b). 

d.  Coordinate  card. 

The  layout  of  the  coordinate  cards  is  the  same  as  for  the  standard  devia¬ 
tion  cards.  Each  coordinate  card  contains  the  fractional  coordinates  X,  Y,  and 
Z  of  an  atom  whose  standard  deviations  are  given  in  the  corresponding  card  of 
the  standard  deviation  deck. 


11.3  Data  Deck 

a.  preliminary  card. 

b.  unit  cell  card. 

c.  all  standard  deviation  cards. 

d.  all  coordinate  cards. 


11.4  CXitput 

Output  may  be  on  logical  tape  9  for  peripheral  punching  or,  if  sense  switch  1  is  down, 
the  cards  are  punched  on-line. 

Each  matrix  card  is  numbered  sequentially  and  identified  with  the  letters  SD  (format 
8E9.4, 14,  3H  SD). 

Each  atomic  coordinate  card  is  numbered  sequentially  and  identified  with  the  letter 
P  (format  8F9.6,  14,  24  P). 

The  two  decks  are  separated  by  a  blank  card. 
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11.5  Limitations 

Die  program  will  accommodate  no  more  than  64  atoms. 


11.6  Special  Subroutines  Called 
None. 


11.7  %)ecial  Tape  Requirements 

Logical  tape  9  is  used  for  peripheral  punching. 


12.  OTHER  FUNCTIONS  AND  SUBROUTINES 
12.1  Function  ARCSIN(X) 

This  program  uses  Hasting’s  Chebyshev  approximation  (sheet  39  of  Ref.  9). 

If  |X  I  >1,  the  program  calls  a  subroutine  ENDJOB  which  may  be  constructed  by  the 
user  to  suit  his  purposes. 

The  program  was  coded  by  B.  A.  Schoomer.'*' 


12.2  Subroutine  RECIP  (AR,  BR,  CR,  ALR,  BER,  GAR,  AA, 

BB,  CC,  COSAL,  COSBE,  COSGA) 

Given  the  sets  of  three  translational  and  three  angular  parameters  for  a  three- 
dimensional  space  lattice,  this  subroutine  finds  the  parameters  for  the  corresponding 
reciprocal  space  lattice. 


Given  Lattice:  Reciprocal  Lattice; 


AR,  BR,  CR,  ALR,  BER,  GAR  AA,  BB,  CC,  AL,  BE,  GA 

AR,  BR,  and  CR  may  be  given  in  any  convenient  units;  ALR,  BER,  and  GAR,  are  given  in 
degrees. 


'•'Diffraction  Brwch,  Optics  Division,  NKL. 
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12.3  Function  DOTPRD  (U,  V,  W,  X,  Y,  Z,  Bl,  B2,  B3,  COSA,  COSB,  COSC) 

'nils_function  calculates  the  scalar  product  of  two  vectors  U  and  V  where  U-U(U,V,W) 
and  V  =  V(X,Y,Z). 

The  vectors  U  and  ^  are  defined  on  a  skew-vector  basis  Bl,  B2,  B3 

BS 


such  that 

U  =  UBl  +  VB2  +  WB3 
V  =  XBl  +  YB2  +  ZB3 
and 

0  •  V  =  UX(B1)2  +  VY(B2)*  +  1Z(B3)*  +  (VZ  +  YW)  B2  B3(ce>»  A)  + 

(UZ  +  fX)  Bl  B3(cos  B)  +  (UY  +  VX)  Bl  B2(cos  C). 

12.4  Function  SCAFAC  (A,  Al,  B,  Bl,  C,  S) 

This  function  calculates  an  atomic  scattering  factor  using  the  analytic  approximation 
of  Vand,  Eiland,  and  Pepinsky,(6). 

The  scattering  factor  f(s)  is 

f(s)  =  A  exp(-Al  S*)  +  Bexp(-Bl  S*)  +  C. 

The  quantities  A,  Al,  B,  Bl,  and  C  are  taken  from  the  tables  of  Forsyth  and  Wells  (2). 

S  is  (sin  o')/  \. 

12.5  Subroutine  MTXMUL  (L,  M,  N,  A,  B,  C) 

This  subroutine  performs  matrix  multiplication  according  to  the  equation 

AB  =  c 

where 

the  number  of  rows  in  A  is  L 
the  number  of  columns  in  A  is  M 
the  number  of  cqlumns  in  B  is  N. 

(As  listed  in  this  report,  the  dimensions  of  all  matrices  are  3x3;  however,  these  may  be 
changed  to  any  other  values  needed.) 
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12.6  Function  INTG(A)  | 

nils  program  converts  floating  point  integers  to  fixed  point  integers  without  trunca¬ 
tion  error.  This  is  done  in  a  manner  such  that  the  integer  value  zero  is  not  printed  with 
a  leading  minus  sign. 
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FORTRAN  LISTINGS  OF  GENERAL  UTILITY  PROGRAMS  APPLICABLE 
TO  CRYSTAL  STRUCTURE  ANALYSIS 


In  the  form  presented  in  this  appendix,  the  FORTRAN-language  programs  are  writ¬ 
ten  for  use  with  the  IBM  7090  IB  Monitor  system. 


Program  Page 

Least-Squares  Cell  Parameter  33 

Subroutine  DET  34 

Subroutine  PARAM  35 

Subroutine  ERRPRM  36 

Subroutine  ERRREL  37 

Subroutine  OUTPRM  37 

Quasi  Normalization  of  Structure  Factors  39 

Rational  Dependence  41 

Subroutine  OUTRAT  43 

Sigma-2  Listings  44 

Subroutine  OUTSIG  45 

Subroutine  REJVEC  46 

Subroutine  SIGVEC,  Triclinic  46 

Subroutine  SIGVEC,  Monoclinic  47 

Subroutine  SIGVEC,  Orthorhombic  48 

Triple  Product  Summation,  Centrosymmetric  48 

Triple  Product'  Summation,  Noncentrosymmetric  50 

Structure  Factor  Calculation  52 

Subroutine  OUTSFC  56 

Interatomic  Distance  and  Angle  57 

Subroutine  INTERP  60 

Subroutine  OUTBND  60 

Least-Squares  Line  and  Plane  Fitter  63 

Point  to  Peak  Distance  Calculation  67 

Form  Factors  for  Busing  Least  Squares  68 

Variance -Covariance  and  Atomic  Parameter  Input  for 
Busing  Function  and  Error  69 
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Program  (Cont’d) 

Page 

Function  ARCSIN(X) 

70 

Subroutine  RECIP 

71 

Function  DOTPROD 

71 

Function  SCAFAC 

71 

Subroutine  MTXMUL 

71 

Function  INTG(A) 

71 

Least-Squares  Cell  Parameter 

DIMENSION  NAME(10)»M(3tl50).W(150)»S(150) .S2( 150) .SC( 150 ) .  02691 

1  H(3t3tl50) t6(3«3) tSG(3«3).BE(3»3)fAL(3«3) (BO) »A { 3 ) »SA( 3 ) »SAL ( 3» 302691 

2  )fEB(3)  02691 

ITP»5 

JTP«6 
50  MM«1 

READ  INPUT  TAPE  I  TP , 1000 ,N1 # U 1ST tYAM, ( NAME ( I ) . I «1 » 10 ) 

IF(Nl)500,500tl00 

100  READ  INPUT  TAPE  I  TP «2000 tM ( 1 1 MM ) »M(2 tMM) »M { 3>MM ) •$ ( MM ) tW ( MM) 


IF{XABSF(M(  l.MM) ) +XABSF ( M( 2 .MM) ) +  XABSFtMt  3.MM ) ) ) 140 . 140 » 110  02691 

110  S2(MM)  =  4,0*(S(MM)  **2 )  02691 

IF(W(MM) )120(120«130  02691 

120  W(MM)  ■  SIMM)  02691 

130  MM  *  MM  +1  02691 

GO  TO  100  02691 

140  MM=MM-1  02691 

IF(YAM)145.145.141  02691 

141  YAM  *  YAM**2  02691 

DO  142  I  =1.MM  02691 

142  S2(I)a  S2(1)/YAM  02691 

145  NN  <  1  02691 

DO  147  I»1.3  02691 

DO  147  J-1.3  02691 

DO  147  K-l.MM  02691 

147  H(  I  tJflO^O.O  02691 

GO  TO(310.290.240.,210.180.150)»N1  02691 

150  N  «  1  02691 

DO  160  I  «1.MM  02691 

160  H(l»l»n«  M(lil)**2+  M(2.I)**2+M(3»I)**2  02691 

GO  TO  325  02691 

180  N>2  02691 

NN«2  02691 

DO  190  I  -l.MM  02691 

H(2.2.I)«  M(3.I)**2 

190  H(l»l»n«M(2*I)**2  +  M(1.I)**2  +  M(2»I)*  M(1»I) 

GO  TO  325  02691 

210  N«2  02691 

DO  220  I«liMM  02691 

H(2t2.n«  M(3.I)*«2 
220  H(1»1»I)«  M(2.1)»*2  M(1.I)**2 

GO  TO  325  02691 

240  N  •  3  02691 

GO  TO  315  02691 

290  N  ■  4  02691 

GO  TO  315  02691 

SIO^N  ■  6  02691 
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I 


I 


315  DQ-  320  I  »  1,3  02691 

DO  320  J  >  1,3  02691 

DO  320  K  >  1,MM  02691 

320  ■  M(I,K)*  M(J,K)  02691 

325  DO  326  I  >  1,3  026.91 

DO  326  J  »  1,3  02691 

G(I,J)  -  0,0  02691 

SG( I , J)  >  0,0  02691 

BEI I,J)  -  0,0  02691 

326  AL ( 1 , J)  >  0,0  02691 

IF(MM-N)327,327,328  02691 

327  WRITE  OUTPUT  TAPE  JTP,3000 

60  TO  500  02691 

328  CALL  PARAM(N,H,W,G,SG,F,SC,MM,S2,LL)  02691 

DO  345  I  ■1,LL  02691 

DO  345  J  sl,LL  02691 

IF(I  -J) 340,330,340  02691 

330  Bd)  =  SORTF(G(I,ln  02691 

GO  TO  345  02691 

340  BE(I,J)=  G(  I  ,J)/S0RTF(G(  I,n*G(  J,J)  )  02691 

345  CONTINUE  02691 

IF(N-2)346,347,362 

346  B(2)»B(1) 

G(2,2)»G(1,1) 

SG<2,2)=SG(1,1)  02691 

GO  TO  360  02691 

347  GO  70(360,350) ,NN  02691 


350  BE(1,3)»  0,5 

BE(3,1)>BE(1,3) 
G(l,3)  a  0,5*G(1*1I 
G(3,l)  a  G(l,3) 

360  B(3)aB(l) 

G(3,3)aG(l,l) 


SG(3,3)*SG(1,11  02691 

362  E8(l)=90.0  -57 .29578*ARCSIN ( BEI 2 ,3 ) ) 

365  EBI2)=90.0  -57.29578*ARCSIN ( BE ( 1 ,3 ) ) 

368  EB(3)=90.0  -57.29578*ARCSINI BEI 1 ,2) ) 

369  CALL  PECIP(BI1),6I2),BI3),EBI1),EBI2),EB(3),AI1),AI2),AI3),  02691 

1  ALI2,3)  ,ALI1,3),AL(1,2))  02691 

CALL  ERRRELIG,SG,BE,A,AL,N,SA,SAL,V,SVR,LL)  02691 

VR  a  1.0/V  02691 


CALL  OUTPRMI A,AL,NAME,B,EB,VR,SAtSAL,SVR,LIST,MM,N,LL,F,SC,S2,W,M 


1,N1,JTP) 

GO  TO  50 
500  CALL  EXIT 

1000  F0RMATI2I2,E8,5,10A6)  02691 

2000  FORHATI3I4,2E10,5  )  02691 

3000  FORMATI 1H1,3X,98H  THERE  ARE  NOT  ENOUGH  DEGREES  OF  FREEDOM.  UNIT02691 
1  CELL  PARAMETERS  AND  ERRORS  CANNOT  BE  CALCULATED*)  02691 

END  02691 


Subroutine  DET 


C  SHARE  DISTRIBUTION  NO*  635 

SUBROUTINE  DET I  A, ALPHA, N, BET) 


DIMENSION  AI6,6) 

BETal. 

1F(N-1)400,200,300  02691 

200  BET  a  AI1,1) 

GO  TO  400  02691 

300  A(l,l)-  A(l,l)-  ALPHA  02691 

6  DO  15  Ia2,N  DET  0009 

A(1,I)>AII,I)-XLPHA  DET  0006 

II-l-l  OET  0007 


I 
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7 

DO  15  J-1«II 

DET 

0008 

8 

IF  (A(I»J))9.15i9 

DET 

0009 

9 

IF(ABSF(A( J.J) )-A0SF(At I .!»)» 11*10.10 

DET 

0010 

10 

R>A(I*J)/A(J*J) 

DET 

0011 

60  TO  130 

DET 

0012 

11 

R«A( JtJ)/A( I tJ) 

DET 

0013 

00  12  K>1*N 

DET 

0014 

BoAIJ.K) 

DET 

0015 

A(J*K)-A(1.K) 

DET 

0016 

12 

A( I.K)>B 

DET 

0017 

BET»-BET 

130 

JJ=J+1 

DET 

0019 

13 

DO  14L=JJ*N 

DET 

0020 

14 

A(ItL)>A(I*L)-R*A(J*L) 

DET 

0021 

15 

CONTINUE 

DET 

0022 

16 

DO  20  I-1»N 

DET 

0023 

20 

BET>BET*A(  I  .n 

400 

RETURN 

02691 

END 


Subroutine  PARAM 


SUBROUTINE  PARAM ( N *H*W .GtSOtF *SC *MM*S2*tL > 

02691 

DIMENSION  H(3*3*150).W(150).6(3*3)*SG(3.3)*SC(150)*S2(150) *013* 

3) .02691 

1 

:  A(6*6)*B(6)«AA(6.7)*X(6«1) 

02691 

IF(N  -3)362.365*365 

02691 

362 

LL  «  N 

02691 

GO  TO  330 

02691 

365 

LL  «  3 

02691 

330 

00  360  I  =  l.LL 

02691 

DO  360  J  *  l.LL 

02691 

IF(  I-J)350.340.3;)0 

02691 

340 

0(1.1)  >  1.0 

02691 

GO  TO  360 

02691 

350 

0(1. J)  =  0.0 

02691 

360 

CONTINUE 

02691 

367 

DO  440  I  =1.LL 

02691 

DO  440  J  =I.LL 

DO  440  K  »1.LL 

02691 

DO  440  L  *K.LL 

02691 

IF( I  -  J)380.370.3a0 

02691 

370 

JJ  =  I 

02691 

GO  TO  390 

02691 

380 

JJ  =1+I+J 

390 

IF()C  -  L)410.400.410 

02691 

400 

KK  »  K 

02691 

GO  TO  420 

02691 

410 

KK  «1+K+L 

420 

A( JJ*KK)-0.0 

02691 

B(JJ)  «  0.0 

02691 

DO  430  1  1  =  1. MM 

02691 

A(JJ.KK)  «  A(JJ.KIC)  +  W(II)*  H(I.J.II)»H(IC.L.II  ) 

02691 

430 

B(JJ)  =  B(JJ)  +  W(II)*  S2(II)*  H(I.J.II) 

A(JJ*)(K)  =(2.-D(I.J))  4(2.0  -  D(K*L)  )»A(JJ*K)() 

02691 

B( JJ)  »  B( JJ)4(2.0  -  D(  I.J)  ) 

02691 

AA(JJ*K)()=  A(JJ*KK) 

02691 

440 

AA(JJ.N+1)  »  B(JJ) 

IF(N-1)445.445.446 

02691 

445 

G(1*1)>B(1)  /A(l.l) 

60  TO  480 

446 

CALL  MATS(AA*X*N.l ) 

02691 

DO  470  I  «  l.LL 

02691 

DO  470  J  «  l.LL 

02691 

IF(1  -  J)  460.450.460 

02691 
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450  G(1»I)«  X(Iil)  02691 

GO  TO  470  02691 

460  K  *1+I+J 

G(I>J)  -  X(Ktl)  02691 

470  CONTINUE  02691 

480  CALL  ERRPRM(WtHf F>SC>SG<N«MM*G*S2tA»LLtD)  02691 

RETURN  02691 

END  02691 


Subroutine  ERRPRM 

SUBROUTINE  ERRPRM ( W tH,F.SC»SG*NtMM*GiS2 tAtLL >0 )  02691 

DIMENSION  M(3t3«150)«SC(150)*SG(3»3)»G(3t3)>S2(150)tA(6f6)*AA(6*6  02691 
1  1 tW( 150) •SK(6I t0(3<3)  02691 

F  *  0.0  02691 

DO  490  L  »  1»MM  02691 

SC(L)  »  0.0  02691 

DO  480  1  =  l.LL  02691 

DO  480  J  =  1 ,LL  02691 

480  SC(L)=SC(L)+H( I .J.L)»G( I .J)*(2.0-0( I.J) )  02691 

490  F  =  F  +  W(L)*<SC(L)  -  S2(L))**2  02691 

V  «  MM  -  N  02691 

IF(N.  -1)500.500.505  02691 

500  SG(1.1)»  F/V  02691 

GO  TO  680  02691 

505  00  506  I»1.N 
00  506  J»1.N 

506  AA( I . J)  a  A( I .J) 

CALL  0ET( AA.O.O.N.TED) 

N2  a  N  -1  02691 

DO  640  Kai,N  02691 

DO  570  1*1. N  02691 

DO  570  Jal.N  02691 

1F(IC-1)520.570.510  02691 

510  II  a  I  02691 

GO  TO  530  02691 

520  11=  I  -  1  02691 

530  JF(K  -  J) 550.570.540  02691 

540  JJ  a  J  02691 

GO  TO  560  02691 

550  JJ  a  J  -1  02691 

560AA(II.JJ)  a.A(I.J)  02691 

570  CONTINUE  02691 

CALL  OET( AA.0,0.N2.EDT)  02691 

640  SK()C)  a  (F*E0T)/(TED*V)  02691 

DO  670  I  a  l.LL  02691 

DO  670  J  a  i,ll  02691 

IFII  -  J)660.650.660  02691 

650  SG(I.I)  a  sxd)  02691 

GO  TO  670  02691 

660  K  -l+I+J 

SG(I.J)  a  SK(K)  02691 

670  CONTINUE  02691 

680  RETURN  02691 

END  02691 
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Subroutine  ERRREL 

SUBROUTINE  ERRREL ( GtSGiBE i A* AL»N tSAtSAL tV«SVR*LL )  02691 

DIMENSION  GO«3)  •SGOtS)  tBEOtS)  >A(3)  «AL I  3 •  3  )  •  SA ( 3  )  tSAL (  3 » 3  )  *0 (  3 )  02691 
V  «  6(1*11*  G(2*2)*  G(3.31  +  2.0*  6(2*31*  G(l«31*  6(1*21-6(1*11*  02691 

1(6(2*31  **21-  6(2*21*(G(1*31  **21-  6( 3*3 1* ( 6 ( 1  * 2 1  **21  02691 

V*S0RTF(V1 

SZ  >  0.0  02691 

00  5  1»1*LL 
5  U(I1«0.0 

IF(N-1110*10*20  02691 

10  SA(11*(SG(1*11*(A(11/G(1*11 l**21/2.0  02691 

GO  TO  380  02691 

20  DO  310  I«1*LL  02691 

00  310  J  -1*LL  02691 

DO  310  K  >1*LL  02691 

IF(I  -J1100*310*100  02691 

100  IF(I  -K1200*310*200  02691 

200  1F(J-K1  300*310*310  02691 

300  U(I1  ■  2.0*  SG( J*Kl-( (G(J*K1/(6(J*J1*  6( K.ICl 1 1 **2 1  * ( S6 ( K*K 1  * (6 ( J*02691 

1  J1  **21+  SG(  J*J1*(6(IC*K1  **21  1  02691 

SZ  -  SZ  +(6(1*11*  U(I1 1*(G(1*I1*(G(J*K1  **21+  6 ( I  * J 1*6( I *K 1*  S0RTF02691 
1  (G(J*J1*G(K*K11 1  02691 

310  continue  02691 

SZ  «  SZ/(2.0  *(V**411  02691 

0  «  1*0  02691 

IF(N-41330*320*330  02691 

320  0  >  -0  02691 

330  DO  370  I  «  1*LL  02691 

00  370  J  »  1*LL  02691 

DO  370  K  «  l.LL  02691 

IF(LL-21335*360*335 

335  IF( I-J1340*370*340  02691 

340  IF(  I-IC1350*370*350  02691 

350  IF( J-K1360*370*370  02691 

360  SA(I1  »((  0*(6(J*K1/(A( I1*(V**2I 1 1**21*  U( 1 1  +  (A(I1**21*1  SG(I*I102691 
1/(G(I*I1  **21+  2.0*SZll/2.0 

SAL(J*IL1«  (r.O  -  AL(J*K1**2  1*(6(I*I1*  U(I1  +  (BE(I*J1**2  +( ( BE( I *02691 

1  IC1**21*(1.0-  BE( I*J1**21*(AL(J*K1**211/(1.0  -  BE(I*K1**211  *  G(J*02691 

2  Jl*  U(J1  +  (BE(I*H1**2  +( (BE( I*Jl**21*( 1.0  -  BE ( 1 *IC  1**2 1* ( AL ( J*02691 

3  Kl**211/(1.0  -  BE(I*J1**21  1*  G(K*ICl*  U(  1C  1  1  / ( 2 .0* ( V**2 1  1  02691 

370  CONTINUE  02691 

380  SVRsSZ*2.0  02691 

420  DO  430  I  «  1*3  02691 

430  SVR  »  SVR  +  SG(I*Il/(  6(1*11**21  02691 

DO  435  I  ■1*LL  02691 

435  SAdl  •  S0RTF(SA(I11  02691 

440  SVR*  SQRTF ( SVR/ ( (V**2 1*2.1 1 

RETURN  02691 

END  02691 


Subroutine  OUTPRM 

SUBROufiNE  0UTPRM(A*AL*NAME*B*EB*VR*SA*SAL*SVR*LIST*MM*N*LL*F*SC*  02691 
1S2*W*M*N1*JTP1 

DIMENSION  NAME(101*M(3*1501*W(150l*S2(1501*SC(1501 * AL ( 3 *3 1 *B( 3 1  *  02691 

1A(31*SA(31*SAL(3*31*EB(31*ALD(31 *SAL0(31  02691 

lF{N-4)95  .50*50 
50  DO  100  I  -1*2 

IL  ■  1  +1  02691 

DO  100  J  -IC*3  02691 

L  ■  6-I-J  02691 

60  ALD(L1*  90.0  -  57.295780*  ARCSIN (AL( I tJll 
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70  D  ■  SQRTF(SAL(I«J) )  02691 

100  SALD(L)»  57.295780*ABSF( ABSF(ARCSIN(AL( I tJ)+D) ) -ABSF ( ARCS  I N ( AL ( I . J02691 

1)-D)))  02691 

95  WRITE  OUTPUT  TAPE  JTPf 1000 * INAME 1 1) *  I • 1 . 10 » tMM 

DO  170  1»  l.N  02691 

GO  70(101. 102t  961 .LL 

96  IF(N-4)103i97.103 

97  GO  70(110.120.130.1601.1 

101  GO  70  no  02691 

102  GO  70(110.1251.1  02691 

103  GO  70(110.120.130.140.150.1601.1 

130  WR17E  0U7PUT  7APE  J7P . 1100.B( 3 1 . A ( 3 1 .SA( 3 1 

GO  70  170  02691 

120  WR17E  0U7PU7  7APE  JTP . 1200 .B( 2 1 . A( 2 1 .SA ( 2 1 

GO  70  170  02691 

125  WRI7E  0U7PU7  7APE  JTP . 1100.B( 2 1 . A ( 2 1 . SA ( 2 1 

GO  TO  170  02691 

no  WRITE  OUTPUT  TAPE  JTP . 1300.B ( 1 1 . A( 1 1 . SA ( 1 1 

60  TO  170  02691 

160  WRITE  OUTPUT  TAPE  JTP. 1400.EB( 31 .AUDI  3 1 .SALD( 3 1 

GO  TO  170  02691 

150  WRITE  OUTPUT  TAPE  JTP . 1500 .EB( 21 .ALD( 2 1 .SALD ( 2 1 

GO  TO  170  02691 

140  WRITE  OUTPUT  TAPE  JTP . 1600. EB( 1 1 .ALD( 1 1 .SAUD ( 1 1 

170  CONTINUE  02691 

WRITE  OUTPUT  TAPE  JTP . 1700 .VR.SVR 
WRITE  OUTPUT  TAPE  JTP. 1800 

GO  T0(500,510.520.530.540.550) .N1  02691 

500  WRITE  OUTPUT  TAPE  JTP. 2400 

GO  TO  175  .  02691 

510  WRITE  OUTPUT  TAPE  JTP. 2500 

60  TO  175  02691 

520  WRITE  OUTPUT  TAPE  JTP. 2600 

GO  TO  175  02691 

530  WRITE  OUTPUT  TAPE  JTP. 2700 

GO  TO  175  02691 

540  WRITE  OUTPUT  TAPE  JTP. 2800 

GO  TO  175  02691 

550  WRITE  OUTPUT  TAPE  JTP. 2900 

175  IF(LIST1260.180.260  02691 

180  R  »  MM  -  N  02691 

SO  ■  S0RTF(F/R1  02691 

S02  *  SO*  2.5  02691 

LI  ■  MM/20  02691 

L2  »  MM  -  Ll*20  02691 

L3  =  0  02691 

185  L3  *L3+1  02691 

IF(L1  -  L31200.190.190  02691 

190  K2  •  L3«20  02691 

K1  »  K2  -19  02691 

GO  TO  220  02691 

195  L3*L3-t-l  02691 

200  IF(L2)260.260.210  02691 

210  K1  »  Ll*20  +1  02691 

K2  =  K1  +  L2  -1  02691 

220  NM  =  0  02691 

WRITE  OUTPUT  TAPE  JTP. 2000. L3.S0 

DO  250  1  .  K1.K2  02691 

DEL  «(S2(n-SC(n)*S(3RTF(W(n)  02691 

IF(SD2-ABSF(DEL) 1230.230.240  02691 

230  NM  >  1  02691 

WRITE  OUTPUT  TAPE  JTP .2100.M( 1 . I ) .M( 2. I ) .M ( 3 . I ) .S2 ( I ) .SC( I ) .W( I ) . 
lOEL 

GO  TO  250  02691 
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240  WRITE  OUTPUT  tape  JTP • 2200*M ( 1 » I ) tMI 2 1 1) tMI 3> I  I «S2 ( 1 ) »SC ( I ) iW ( I ) » 
lOEL 

250  CONTINUE 

IF(NM)256*256i255 

255  WRITE  OUTPUT  TAPE  JTPf2300 

256  IF(L1-L3)260il95.185 
250  RETURN 

lOuO  FORMAT! 1H1/1H4»36X.45HUNIT 
lS//29XflOA6//46X I23HNUMBER 
1100  FORMAT! lH0.26Xf5H  C* 

1200  FORMAT! lH0t26X.5H  B* 


02691 
02691 

02691 
02691 

CELL  PARAMETERS  WITH  STANDARD  DEVI AT10N02691 

OF  REFLECTIONS  »iI6//l  02691 

=»F10.7«22Xt4H  C  *.F9.5»6H  +OR-#F8.5/)  02691 

=»F10.7*22Xf4H  B  »»F9.5.6H  +OR-.F8.5/)  02691 


1300 

1400 

1500 

1600 

1700 

1800 

2100 

2200 

2000 


FORMAT! 1H0.26X»5H  A*  = tF10.7«22X*4H  A  =.F9.5.6H  +OR-.F8.5/)  02691 

FORMAT! lH0t23X»8HGAMMA*  = *F10.5 . 19X .7HGAMMA  =.F9.4.6H  +OR-»F8 tA/ ) 02691 

FORMAT! 1H0.23X.8H  BETA*  » iFlO.5* 19X»7H  BETA  ».F9.4t6H 
FORMAT! lH0t23X.8HALPHA»  * .FIO. 5» 19X .7HALPHA  =»F9«4»6H 
FORMAT! 1H0.64X«3HV  =.F11,4,6H  +OR-,F10.4//// ) 

FORMAT!33X.48HUNITS  ARE  1 /! ANGSTROMS  I .  ANGSTROMS*  AND 
FORMAT! lH0«llX«3l6*6X«E14.7*5X*E14*7t3X*E12.5>6X*E12.5*3H 
FORMAT! lH0*llX«3I6*6XtE14.7«5X.E14*T*3X*E12.5>6X*E12*5 I 


+OR-.F8. 4/ 102691 
+OR-*F8.4/ >02691 
02691 

DEGREES)  02691 
*)  02691 

02691 


F0RMAT!1H1.45X.25H0BSERVE0  AND  CALCULATED  S.25X.5HPAGE  . 1 3///27X. 10269 1 
18HS  =  !2*SIN! )/L)**2.10X.21HSTAN0ARD  DEVIATION  *  .E 12.5 ///17X . 13HH0269 1 
2  X  L*10X.6HS!0BS)<13X>6HS!CAL).12X*1HW.  9X.16HIDELTA  S)*!W02691 


3**l/2)) 

2300  FORMAT! 1H0/9X.58H* 
lEVIATION) ) 

2400  FORMAT! 1H0.45X.27HTHE 
2500  FORMAT! 1H0.45X.27HTHE 
2600  FORMAT! 1H0.45X.27HTHE 
2700  FORMAT! 1H0.45X.27HTHE 
2800  FORMAT! 1H0.45X.27HTHE 
2900  FORMAT! 1H0.45X.27HTHE 
END 


THIS  DEVIATION  IS  GREATER  THAN  2.5* ! STANDARD 


CRYSTAL 

CRYSTAL 

CRYSTAL 

CRYSTAL 

CRYSTAL 

CRVSTAL 


IS 

IS 

IS 

IS 

IS 

IS 


TRICLINIC 

MONOCLlNIC 

ORTHORHOMBIC 

tetragonal 

HEXAGONAL 

CUBIC 


02691 

002691 

02691 

02591 

02691 

02691 

02691 

02691 

026911 

02691 


Quasi  Normalization  of  Structure  Factora 

CSFQNOR  QUASI  NORM.  OF  STRUCT.  FACTORS  USING  F  TAPE  INPUT 

DIMENSION  NME!12).  AS ! 20 ) . ASl 1 20 ) ,BS! 20 ) .BSl ! 20 ) . CS ! 20 > . AN ! 20 ) .  02691 


1  HI! 500) .H2! 500) .H3 ! 500 ) .Ml ! 500 ) 
ITP»5 

JTAPE  «  15 

ITAPE  =  16 

REWIND  ITAPE 
REWIND  JTAPE 
READ  INPUT  TAPE  ITP. 

1  lOOO.NOIB.NASP.NF.BTl.EX 

READ  INPUT  TAPE  ITP. 

1  40.!AS!I).AS1!I).BS!I).BS1! 

IF!NF)5.5.1 

1  DO  2  1*1. NF 

2  CALL  FLSXPD! ITAPE) 

5  IF!LIB)10»30.10 

10  READ  INPUT  TAPE  ITAPE. 1200. !NME! 
30  NO  ■  0 

N1  «  NOIB/500 

WRITE  OUTPUT  TAPE  JTAPE. 740. NOIB 
100  IF!N1-N0)400. 150.200 
150  N4  ■  NOIB  -  Nl*  500 
IF!N4)400 *400.210 
200  N4  >  500 


M2!500) .M3! 500) .AVE0I500) .SI ! 500)02691 


02691 

02691 

02591 

02691 

.LIB 

)»CS!I).AN!I).I>1.NASP) 


) .1*1.12) 

02691 

02691 

02691 

02691 

02691 

02691 

02691 
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210  READ  INPUT  TAPE  ITAPEi 1100 > (HI ( I ) tH2 I  I ) tHS ( 11 • AVEQ ( 1) »SI ( 1) t !■  It  02691 
1  N4)  02691 

DO  300  I  ■  ltN4  02691 

M1(I)»  INTG(Hl(n)  02691 

M2II)>  INTG(H2(I))  02691 

M3( I INTG(H3(I ) )  02691 

SI ( I )>  SI ( I )••2  02691 

SIG  -  0.0  02691 

DO  250  J  »  ItNASP  02691 

2  50  SIG»SlCi+AN(  J)»(  AS(  J)*EXPF<-AS1(  J)*SI(  1  )  )+  BS  (  J I  *EXPF  ( -BSK  J )  *SI  ( 1)02691 

1)+CS(J))**2  02691 

300  AVEQ{  I)  =  (AVE0(  n*  EXPF ( BT  1*( SI  ( 1 ) **EX )  ) /S IG ) -  1.0  02691 

WRITE  OUTPUT  TAPE  JTAPE .680* (Ml ( I ) tM2 ( 1) tM3 ( I ) t AVEO ( I ) t I  >ltN4)  02691 
NO  »  NO  +1  02691 

GO  TO  100  02691 

400  REWIND  ITAPE  02691 

END  FILE  JTAPE  02691 

REWIND  JTAPE  02691 

PRINT  2000  02691 

CALL  EXIT 

1000  FORMAT(3I10,2E10.4tI10) 

40  FORMAT(5F8.4tF4«0)  02691 

1100  F0RMAT(4F9.2tl8X.F9.6)  02691 

680  FORMAT(3I4tF20.5)  02691 

2000  FORMAT!  16H  JOB  FINSHED, /////////// ) 

740  FORMATd?)  02691 

1200  FORMAT(12A6) 

END 
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Rational  Dependence 

CRATDEP  RATIOiNAL  .IKPE.NJA'K  E  (AFTrR  ULOCK-YANNON 1) 

DIMENSION  NAME(12).  NF  (  2  =;  )  .  KA  (  ?fl  )  ,KH  (  28  )  t 
1ICC(26  )  fEC(29) .MK( 5000) .ML ( 5000) .MH) 5000) .ESOM( 5000 ) .MR ( 5000) 
COMMON  MH.MK.ML.ESOM.MR.KA.KG.KC.jrjF.  lAO.M.IT.NAME.N.J.K.MA.JTP. 

iavesog.jput.nctr.no 


KA(1) 

=  0 

06800 

KA(2) 

=  1 

06800 

KAO) 

=-l 

06800 

KA(4) 

=  2 

06800 

KA(5) 

=-2 

06BO0 

KA(6) 

=  3 

06800 

KA(7) 

=-3 

06300 

KA(8) 

=  4 

06800 

KA(9) 

=-4 

06800 

KA(IO) 

*  5 

06800 

KA( 11  ) 

*-5 

06800 

KA( 12  ) 

=  6 

06800 

KA( 13  ) 

=  -6 

06800 

KA( 14) 

=  7 

06800 

KA( 15  ) 

=-7 

06800 

KA(16) 

=  8 

06800 

KA( 17) 

=  -8 

06800 

KA( 18  ) 

=  9 

06800 

KA( 19) 

=-9 

06800 

KA(20) 

=  10 

06800 

KA(21  ) 

=-10 

06800 

KA(22  ) 

»  11 

06800 

KA(23  ) 

S-11 

06800 

KA(24) 

=  12 

06800 

KA(25) 

=  -12 

06800 

KA(26) 

=  13 

06800 

KA(27) 

=  -13 

06800 

KA(28) 

=  14 

06800 

ITP  =  5 

JTP  =  6 

KTP=12 

REWIND 

KTP 

50  READ  INPUT  TAPE  I  TP . 8 .  (  N AME ( I ) . I =1 . 12 ) 

READ  INPUT  TAPE  I  TP .4 .NA .N I  flT .NF IN »MI NN .AN .SDMI N» JPUT 

60  READ  INPUT  TAPE  KTP.5.NO 

61  READ  INPUT  TAPE  K.TP  .6  .  (  MH  (  J  )  .MK  (  J )  .ML  ( J  )  .ESQMI  J  )  .  J=  1  .NO  ) 

REWIND  KTP 

SUMESO=0.0 

DO  62  I=liNO 

ESOMI I )=ESOM( I )+1.0 

62  SUMESQ=SUMESQ+ESQM( I  ) 

BNO=NO 

AVESO=SUMESC!/BNO 
70  MA*NINT 
NCTR=-1 

WRITE  OUTPUT  TAPE  JTP . 15 . ( NAME ( J ) » J  =  1 » 12 ) .N I  NT .NF IN . AVESO. AN.M I NN  . 


INO.SDMIN 

80  1F(NA)83.90.90  06800 

83  AN»1.414214*AN  06800 

90  DO  95  J=  1.28 

92  KB(J)alCA(J)  06800 

95  <C(J)=KA(J)  06800 

180  DO  440  N«1.MA  06800 

190  DO  440  J«1.MA  06800 

200  DO  440  IC«1.IMA  06800 

N«N 

JmJ 

K*K 


IF(KA(N))  440.2000.202 
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2000 

IF(KB( J) )440i2010.202 

2010 

IF(KC(<) )44Ci440.202 

202 

MAD=MA+1 

DO  207  L0=ltMA0 

NF(LD)=0 

207 

EC(LD)=0. 

2ino 

DO  2400  L*ltN0 

KSUM»KA ( N ) *MH ( L ) +KB ( J ) •MK ( L ) +KC ( < ) *ML ( L ) 
MT=XMODF(<SUM»MA) 

IF(MT)2200.2300i2300 

2200 

MT=MT+MA 

2300 

MT=MT+1 

MR(L)=MT 

EC(MT)=ESQM(L)+EC<MT) 

2400 

NF(MT)=NF(MT)-fl 

300 

PMA*MA 

301 

PMA=PMA/2. 

302 

man'=pma 

303 

MAN*MAN+1 

3  04 

IT=MA+2 

310 

DO  430  M=liMAN 

311 

IT=IT-1 

1315 

IF(IT-M)  1320.1321*1325 

1320 

STOP  0066 

1321 

0NF=NF(M) 

JNF=NF(M) 

1322 

BEC=EC(MI 

1323 

GO  TO  1350 

1325 

JNF=NF(M)+NF(IT) 

1330 

0NF=JNF 

1340 

BEC*EC(M)+EC(IT) 

1350 

IF(JNF-M1NN  >430.430. 1355 

1355 

AVESOG=BEC/DNF 

1360 

CHEK=AN/SORTF(DNF| 

IF (CHEK-SDMIN) 1365.1370.1370 

1365 

CHEK*SDMIN 

1370 

PHEK=AVESQ+CHEK 

1360 

PMHEK=AVESO-CHEK 

360 

IF(AVESOG-PHEK)  370.380.  380 

06600 

370 

IF(AVESQG-PMHEK)380.  380.  430 

380 

IF  (MA-3)400. 400. 3193 

3193 

IR  =  0 

06800 

3194 

CAT=2.0 

06800 

3200 

XMA=MA 

06800 

3210 

XLA=XMA/CAT 

O68O0 

3215 

IZ=«0 

06800 

3220 

MX=XLA 

06600 

3230 

YMX=MX 

06800 

3240 

IF  (YMX-XLA)  3350.3250.3350 

06800 

3250 

XKE=XA(N) 

06800 

3260 

AAD=XKE/CAT 

06600 

3270 

MAC=AA0 

3280 

AB0=MAC 

3290 

IF(ABD-AAD)  3350.3300.3350 

06600 

3300 

IZ»I2+1 

06800 

3310 

GO  TO  (3320. 3330. 3340. 430). IZ 

3320 

XICE=KB(  J) 

06800 

3325 

GO  TO  3260 

06800 

3330 

XKE=KC(K) 

06800 

3335 

GO  TO  3260 

06800 

3340 

XKE  «M-1 

3345 

GO  TO  3260 

06800 

3350 

IR-IR+1 

06800 

3355 

IF(MA-IR-3>  400.400.3360 

06800 

3360 

CAT»CAT+1. 

06800 
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3365  IF(CAT-14,0|  3210.3210.400 
4C0  CALL  OUTRAT 
430  CONTINUE 

IFISENSE  SWITCH  2)510.440 

440  CONTINUE  06800 

450  MAxMA+l  06800 


460  IF(NFIN-MA)  470.180.180 

470  READ  INPUT  TAPE  I  TP .4 .NA .N I  NT .NF IN .MINN .AN .SDMI N . JPUT 
IF(NINT)500.500.70 
500  CALL  EXIT 

510  WRITE  OUTPUT  TAPE  JTP . 9 . { NAME ( L2 )  .L2=  1 . 12 ) 

WRITE  OUTPUT  TAPE  JTP . 7 .KA ( N ) .KB ( J ) .KC( K ) .MA 
WRITE  OUTPUT  TAPE  JTP. 16 
GO  TO  500 

8  FORMAT  (12A6) 

4  FORMAT  (4Iie.2E10.4.I10) 

5  FORMAT  (17) 

6  FORMAT(3I4.F20.5) 

7  FORMAT(6HO  A  >I3.6X.3HB  3l3.6X.3HC  3 1 3 ,6X.6HMODE  3I6) 

9  FORMAT  ( 1H1.19X.12A6) 

15  FORMAT(1H2.19X.12A6////20X.16HINITAL  MODULAS  = I  6  .  lOX  .  15HF I NAL  MODU 

ILAS  3t6///2CX.19HAVERAGE  E**2  »F9.5///20X.  13HDEVI AT  IONS  OFF7 

2.3.  66H  STANDARD  DEVIATIONS  FROM  AVE  E*»2  ARE  CONSIDERED  SIGN 
3IFICANT///2OX.  40HONLY  SUBGROUPS  CONTAINING  NO  FEWER  THAN  16, 

4  39H  REFLECTIONS  ARE  CONSIDERED  SI6NIFICANT///20X.  lOHTHERE  ARE  17 
5.  30H  REFLECTIONS  ON  THE  INPUT  TAPE///20X,51HTHE  MINIMUM  ACCEPTED 
6DEVIATION  FROM  AVE.  E**2  IS  F5.3> 

16  FORMAT(1H0.19X.47HCALCULATION  TERMINATED  BY  SENSE  SWITCH  CONTROL.) 
END 


Subroutine  OUTRAT 

DIMENSION  NAME(12).  NF ( 29 ) .KA ( 28 ) .KB ( 28 ) . 

1KC( 28),EC(29) ,MK( 5000) .ML ( 5000 ) ,MH( 5000 ) .ESQM ( 5000 ) .MR ( 5000 ) . 
2LH(40).LK(40) .LL(40).ESOL(40) 

COMMON  MH.MK.ML.ESOM.MR.KA.KB.KC.JNF.MAD.M.IT.NAME.N.J.K.MA.JTP. 

iavesqg.jput.nctr.no 

KET3M-1 

IF( JPUT)75,50,75 
50  NCTR-NCTR+l 

IF(NCTR-(NCTR/9)*9)60,55,60 
55  WRITE  OUTPUT  TAPE  JTP . 9. ( NAME ( L2 ) .L23 1 . 12 ) 

60  WRITE  OUTPUT  TAPE  JTP. 8 .KA ( N ) ,K8 ( J ) »KC( K ) .MA.KET .  AVESOG.JNF 
GO  TO  350 
75  NP*0 
L4=0 

DO  220  L93l,NO 
L9«L9 

IF(MR(L9)-M)200,210,200 
200  IF(MR(L9)-IT)220.210.220 
210  L4«L4+1 

LH(L4)3MH(L9) 

LK(L4)3MK(L9) 

LL(L«)3ML(L9) 

ESQL(L4)«ES0M(L9) 

IF(40-L4)230.230,220 
220  CONTINUE 
230  Jl«L4/40 

J2»L4-J1*40 

J3»J2-20 

IF  ( J3)241.241.243 
241  J4-J2 

GO  TO  245 
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243  J4»J2-2*J3 
GO  TO  247 

245  IF  ( Jl)3O0i3O0»246 

246  J5»J1 

GO  TO  250 

247  J5-J1+1 

250  DO  290  L1=1»J5 
NP-NP+1 

WRITE  OUTPUT  TAPE  JTP . 9 f ( NAME ( L2 ) »L2* 1  *  12 ) 

WRITE  OUTPUT  TAPE  JTP ♦ 7 ,KA ( N ) .KB ( J ) .KC ( K ) *MA *KET »NP , AVESOG * JNF 
L5«L1*40-39 
IF  ( J1-L1)260.270.270 
260  L6=L5+J3-1 
GO  TO  280 
270  L6*L5+19 
280  DO  290  L2=L5.L6 
AO»ESOL(L2)-1.0 
A0A*ESQL(L2)-AVES0G 
A01=ESQL(L2+20)-1.0 
AQA1=ESOL(L2+20) -AVESOG 
WRITE  OUTPUT  TAPE  JTP.IO 

290  WRITE  OUTPUT  TAPE  JTP.  1 1 .LH  (  L2 ) . LK ( L2 ) »LL t L2 ) » AO .AOA.LH ( L2+20 ) . 

lLK(L2+20) .LL(L2+201 .AQl.AOAl 
300  IF  ( J4)360.360.310 
310  IF  (03)320.320.315 
315  L5  =  L6  +  1- 
GO  TO  335 
320  NP=NP+1 

WRITE  OUTPUT  TAPE  JTP . 9 . ( NAME ( L2 ) .U2* 1 . 12 ) 

WRITE  OUTPUT  TAPE  JTP . 7 .KA ( N ) .K8 ( J ) .KC( K ) .MA.KET .NP .AVESOG . JNF 
330  L5*L6+21 
335  L6*L5-fJ4-l 

DO  340  L2=L5.L6 
AO=ESQL( L2)-1.0 
AQA=ESQL(L2) -AVESOG 

340  WRITE  OUTPUT  TAPE  JTP . 12 .LH ( L2 ) . LK( L2 ) .LL ( L2 ) . AO.AQA 
350  RETURN 
360  L4*0 

IF(NO-L9)350.350.220 

7  FORMAT(6HO  A  =I3.6X.3H6  =I3.6X.3HC  » I  3 .6X.6HM0DE  •I6.8X. IIHREMAIN 

IDER  =I6.24X.4HPAGEI5///22X.10HAVE  E**2  =F 10. 5 . 8X  .  IIHNO.  REFL.  =16/ 
2///2(60H  H  K  L  E**2  -  1  E**2  -  A 

3  )  I 

8  FORMAT(6HO  A  =I3.6X.3HB  =I3»6X.3HC  = I  3 »6X .6HM0DE  *  1 6 .8X » IIHREMAI N 

IDER  =I6.24X.6H  ///22X,10HAVE  E*»2  =F 10. 5 .8X . 1 IHNO.  REFL.  =16) 

9  FORMAT  ( 1H1.19X.12A6) 

10  FORMAT  (IH  ) 

11  FORMAT  (2(1X.3I4.7X.F9.5.7X.F9.5.15X) ) 

12  FORMAT  ( 1H0.314.7X.F9.5.7X.F9.5) 

END 


Slgma-2  Listings 


CS16EXC  SIGMA-2  INTERACTION  EXC 

DIMENSION  Ml(2400) .Mai  2400). M3( 2400). AV( 2400) •Kl(2AO0).K2(2§0O)iK3 
1( 2400) .OAVI 2400) 

COMMON  Ml. M2. M3. AV 

ITP«5 

IP  «X1 

L»0 

READ  INPUT  TAPE  ITP. 

1  1002.Z.I2.NOIB 


IF(IZ)  2.  1.2 


02B91 


NAVAL  RESEARCH  LABORATORY 


45 


1  MM»2 
REWIND  IP 

READ  INPUT  TAPEIP.lOOliNRT 
GO  TO  3 

2  MM«1 

READ  INPUT  TAPE  ITPilOOl.NRT 

3  DO  25  I»ltNRT 
20  GO  TO(5t6) »MM 

5  READ  INPUT  TAPE  ITP, 

1  1000tNlfN2iN3tO 

GO  TO  7 

6  READ  INPUT  TAPEI P . 1000 .N1 »N2 tN3»0 

7  CALL  REJVEC(N1.N2,N3*  0) 
IF(Q-Z)25. 10*10 

10  L»L+1 
MltL)»Nl 
M2(L)«N2 
M3(L)«N3 

AV(L)»SORTF(1.0+0) 

25  CONTINUE 
30  DO  160  J«1«L 
40  DO  70  IC«1*L 
50  IFIAVdC)  170*60.60 
60  BAV«AV(K) 

GO  TO  80 
70  CONTINUE 
80  00  110  1«K*L 
90  IF(  AV(n-BAV)110*100*100 
100  BAV>AVI I  ) 

Il-I 

11-0  CONTINUE 

K1(J)«M1(I1) 

K2( J»«N2( m 

K3( J)»M3( Ill 
OAV( J)»AV(I1) 

AV(  Ill»-1.0 
160  CONTINUE 

IF(NOIB) 165*180*170 
165  NOIB>L 

170  CALL  SIGVEC(NOIB*Kl*K2*K3*OAV*L) 
180  CALL  EXIT 

1000  FORMAT(3I4*F20.5I 

1001  FORMATIIT) 

1002  FORMATIEIO. 4*2110) 

EMO 


02691 

02691 

02691 

02691 


02691 


02691 


02691 


02693 


Subroutine  OUTSIG 


SUBROUTINE  0UTSI6(M1*M2*N3*0*N1*N2*N3*  02693 

101*00*N*K1*K2*K3*I*0AVI  02693 

DIMENSION  K1(2400)*X2(2400)*K3(2400)*  02693 

lOAVI 2400) *M1 ( 1000 ) *M2( 1000 ) *M3( 1000) *0( 1000 ) *  02693 

2N1(1000) *N2I 1000) *N3( 1000) *01 (1000) *001 1000)  02693 

COMMON  M1*M2*M3*0*N1*N2*N3*OI*00  02693 

JTP«6 

IF(N)25*25*50  02693 

25  X»1 

WRITE  OUTPUT  TAPE  JTP* 

1  1000*K*K1( I)*K2( I)*K3(I)*OAV(n 

GO  TO  130  02693 

50  J-N/20  02693 

IF(J)120*120*100  02693 
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02693 

02693 

02693 


100  DO  110  Il-l.J 
J3»20*I1 
J2-J3-19 

WRITE  OUTPUT  TAPE  JTP. 

1  lOOOtlltKK  n.K2(n«K3(I)*QAVtn 

110  WRITE  OUTPUT  TAPE  JTP. 

1  2000.IM1(J1).M2(J1)»M3(J1) .QIJl). 

2N1( J1).N2( Jl) iN3( Jl) .OK Jl) »Q0( J1)«J1-J2.J3) 
IF{N-20*J)130.130.120 
120  I2»20*J+1 
I4-J+1 


02693 

02693 

02693 


WRITE  OUTPUT  TAPE  JTP. 

1  1000.I4.Kl(I).K2(II.K3(n.OAV(I  ) 

WRITE  OUTPUT  TAPE  JTP. 

1  2000.(M1IJ1).M2( Jl) .M3(J1) .0(J1). 

2N1( Jl) .N2( J1).N3( Jl) .OK Jl) .OOI J1).J1-I2.N) 

130  RETURN  02693 

1000  F0RMAT(1H1.17X.9HH  <  L .9X.1HE.13X .5HPA6E  . I 2///16X.3I4.4X.F8.5 
1////70H  H{I)  Kd)  Ld)  Ed)  HIJ)  K(J)  L(J)  E(J) 

2Ed)*E(J)*E) 

2000  FORMAT! IHO. I4.2I5.4X.F8.5.1X.3I5.4X.Fe.3.4X.F10*5) 

END  02693 


Subroutine  REJVEC 


SUBROUTINE  REJVEC (N1 .N2.N3 .0) 

C  DUMMY 
100  RETURN 
END 

SUBROUTINE  REJVEC ( N1 .N2.N3 .0 ) 

C  REJECTS  ALL  REFL.  FOR  WHICH  ANY  OF  H.K.L  EQUAL  ZERO 
IF(N1)10.30.10 
10  IF(N2)20.30.20 
20  IF(N3)40.30.40 
30  0—10. 

40  RETURN 
END 


Subroutine  SIGVEC,  Triclinic 


C  TRICLINIC 

SUBROUTINE  SIOVEC (NOIB .K1 .K2 .K3.0AV.L )  02693 

DIMENSION  K1(2400).X2(2400).K3(2400)>  02693 

1QAV(24.00)  .MldUOO)  .M2(  1000  )  .M3 ( 1000)  «0(  1000)  .  02693 

ZNKIOOO)  .N2(  1000)  .NS!  1000  )  .01  (  lOoO)  .QQ(  1000)  02693 

COMMON  M1.M2.M3.0.N1.N2.N3.01.0Q  02693 

DO  170  I=1»N0IB  02693 

N«0  02693 

DO  160  J>1.L  02693 

DO  160  IC»J.L  02693 


IF(K1(J)+IC1(K)-Kld))  10.  40.  10 
10  IF(Kl(J)+Kl(IC)+Kld))  20.  50.  20 
20  IF(Xl(J)-Kl(K)-tCld))  30.  60.  30 
30  IF(ICl(J)-Kl(K)-flCld))160.  70.160 
40  IF(X2(J)+IC2IK)-X2d))160.  80.160 
50  IF(K2(J)+X2(K)+<2d))160.  90.160 
60  IF(K2(J)-)C2IIC)-K2d  ))  160.100.160 
70  IF(K2(J)-K2(K)-H(2d  ))  160.110.160 
80  IF(K3( J)4X3(K)-K3d )) 160.150.160 
00  1F(K3!  J)+)C3!K)*<3)  I) )  160.150.160 
100  IF(K3!J)-IC3(X)-IC3d))  160.150.160 
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110  IF(K3(J)-K3(K)-fK3(I) )  160 #150.160 
150  N»N+1 

M1(N)>K1(J) 

M2(N)>K2( J) 

M3(N)>K3( J) 

0(N)  >  OAV(J)  02693 

N1(N)>K1(K) 

N2(N)<K2(<) 

N3(N)>K3(K) 

Q1(N»=  QAV(K)  02693 

00(N)«  0AV(J)«  QAV(K)»QAV(n 

160  CONTINUE  02693 

17U  CALL  0UTSIGIM1.M2.M3.0.N1.N2.N3. 01*00. N. 

1  K1.K2.K3.I.OAV) 

RETURN  02693 

END  02693 


Subroutine  SIGVEC,  Monoclinic 


C  HONOCLINIC*  SECOND  SETTING 

SUBROUTINE  SIGVEC ( NOI B *K1 *K2 .K3.0AV.L) 
DIMENSION  K1 ( 2400 ) .<2(2400) *K3( 2400) * 
10AV(2400) .M1(1000).M2( 1000 ) *M3 ( 1000 ) .0( 1000) . 
2N1(1000) .N2<1000),N3(  1000 ) .01 ( 1000 ) .00( 1000) 
COMMON  M1.M2.M3.Q.N1.N2.N3.01.0Q 
00  170  I-l.NOIB 


N»0 

DO 

160  J»1.L 

00 

160  <«J,L 

IF( 

XABSF(<2( J) 

♦K2 

(K) 

10 

IF( 

XABSF<K2(J) 

-<2 

(K) 

20 

IF( 

Kl(J)-*-)LI(<) 

-Kl 

(I) 

30 

IF( 

<l(J)+<l(lt) 

+<1 

(1) 

40 

IF( 

<1(J)-K1(K) 

-Kl 

(I) 

50 

IF( 

m(J)-Xl(K) 

+K1 

(I) 

60 

M«1 

IF( 

<3IJ)4<3(<) 

-K3 

(1) 

70 

M«2 

IF( 

K3(J)+K3(m4IC3 

(I) 

80 

IF) 

K3(J)-K3(<) 

-K3 

( I ) 

90 

IF( 

F.3(  J)-<3(IC)+K3 

(I) 

100 

IF( 

K1(JH120,1 

10* 

120 

no 

GO 

T0(140,145) 

.M 

120 

IF( 

ICKK)  1160,130. 

160 

130 

GO 

TO(14S.140) 

,M 

140 

IF( 

<3(J)-K3(K) 

♦K3 

(I) 

145 

IF( 

<3(J)-K3(K) 

-K3 

(I) 

146 

IF( 

Kl(ini60,90.160 

150 

N«N+1 

-XABSF(<2(I ) ) )  10.  20.  10 
-XABSF(<2(I ) ) )160.  20.160 
30*  60.30 
40.  70.40 
50.  80.50 
160.  90.160 

100.150.100 

100.150.100 
146,150.146 
160.150.160 


160.150.160 

160.150.160 


M1(N1*IC1(J) 

M2(N)«K2(J) 

M3(N)-K3(J) 

0(N)  -  QAV(J) 

N1(N)>K1(K) 

N2(N)>K2(K) 

N3(N)-<3<)() 

01(N)«  QAV(K) 

00(N)«  QAV(J)*  QAV(K)»QAV( I ) 

160  CONTINUE 

170  CALL  0UTSIG(M1.M2.M3»0.N1,N2.N3,Q1,0Q*N., 
1  Kl.K2.K3,I,OAV) 

RETURN 

END 


02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 


02693 


02693 

02693 

02693 

02693 
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Subroutine  SIGVEC,  Orthorhombic 


ORTHORHOMBIC 

SUBROUTINE  SIGVEC ( NOI B iKl tK2 •K3*OAV*L )  02693 

DIMENSION  Kl(2400|tK2(2400)tK3(2400)«  02693 

10AV(2400) »M1(1000) ,M2( 1000 ) .M3( 1000 J ♦01 1000)*  02693 

2N1(1000) •N2( 1000) .N3( 1000) *011 1000) tOOI 1000)  02693 

COMMON  M1*M2*M3*0*N1*N2*N3*Q1*00  02693 

DO  170  I»l»NOIB  02693 

N=0  02693 

DO  160  J»1,L  02693 

DO  160  K«J.L  02693 


IF(XA8SF(K1(  J)+K1  (K)  )-XA8SF(ICl(l  )  ) )  100 , 1 10 *  100 
100  IF(XABSF(K1(  J)-K1(X))-XABSF(ICUI  )  ) )  160. 1 10 *  160 
110  IF(XABSF(K2( J)+K2(K) )-XABSF<IC2II ) ) ) 120. 130 . 120 
120  IFIXABSF(K2( J)-K2(X) )-XABSF(K2(l ))) 160. 130. 160 
130  1F(XABSF(K3( J)+K3(K) )-XABSF(K3(I ))) 140. 150. 140 
140  I F ( XABSF (K3(J)-K3(K)) -XABSF ( K3 ( 1  ) ) ) 160 *  1 50 . 160 


150  N*N+1 

M1(N)=K1(J) 

M2(N)=K2(J) 

M3(N)=K3(J) 

OIN)  =  OAV(J)  02693 

Nl(N)=Ki(K) 

N2(N)»K2(K) 

N3(N)=K3(K) 

01(N)»  QAV(K)  02693 

COIN)*  OAV(J)*  0AV(K)*0AV( I  ) 

160  CONTINUE  02693 

170  CALL  0UTSIG(M1.M2.M3.0.N1.N2.N3.01.00.N. 

1  K1.K2.K3.I .OAV) 

RETURN  02693 

END  02693 


Triple  Product  Summation,  Centrosymmetric 

CORTHTP  OPTH'^RHOMBIC  TRIP.  PROD.  PROG.  FOR  CFNTRO.  CRYST. 
DIMENSION  E(30.10.30) 

ITP  =  5 
JTP»6 
)CTP»9 

REWIND  KTP 
B»1.0E-21 

READ  INPPT  TAPE-  I  TP.  1000  .NA.  IM.  JM.KM 
1000  FORMAT(6I10) 

DO  75  1=1. IM 
DO  75  J-l.JM 
DO  75  K-l.KM 
75  E(  I*J*K)=B 

READ  INPUT  TAPE  KTP.20Q0.NOIB 
2000  FORMAT! 17) 

DO  8  NO-l.NOIB 

READ  INPUT  TAPE  ICTP.2 100. 1  » J.K.A 
2100  FORMAT{3I4.F20.5) 

IF(I)1.8.2 

1  I=-I 

2  1F(J)3*8*4 

3  J"-J 

4  IF(K)5.8*6 

5  K—K 

6  E(I.J.X)=A 
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8  CONTINUE 
REWIND  KTP 
AN»NA 

SF«SORTF(AN**3)/8.0 

JJJ»0 


C 

9 


MN«0 

START  NEW  TP  CALC. 

read  input  TAPE  I  TP » 1000 «M1 »M2 .M3 »N1 *N2 »N3 

I1«XABSF(M11 

I2=XABSF(M2) 


I3*XABSF(M3) 

J1=XABSF(N11 


J2=XABSF(N2) 

J3=XABSF(N3) 

L1»M1+N1 

L2«M2+N2 


L3*M3+N3 

K1=XABSF(L11 

K2=XABSF(L2) 

K3=XABSF(L3) 

IF(I1  +  I2+I3+'J1+J2+J3)  10.10.50 

10  MN»MN+1 
EEE*1.0 
COR=0.0 

GO  T0(55.630) .MN 
50  IF(Ml)11.2t.ll 

11  IF(M2)12.20.12 

12  IF(M3  )  13.20.13 

13  1F(N1I1A.20.14 

14  IF(N2)15.20.15 

15  IF(N3)17.20.17 

17  IF(L1)18. 20.18 

18  IF(L2)19.20.19 

19  IF(L3)25.20.25 

20  EEE=1.0 
C0R=0.0 


25  EEE»SORTF(  (Ed  1 . 1 2 . 1 3  )  •*■1 .0  I  •  (  E  (  J1 .  J2  .  J3 )  *1.0 1  *  (  E  ( K1  .K2  .K3 ) +  1.0 )  ) 
C0R=  (E(I1.I2.I3)+E!-J1.J2.J3)+E(K1.IC2.K.3»+1.0)  /S0KTF(  AN  ) 

C  CALCULATE  TRIPLE  AVERAGES 
55  RN=0.0 
TP»0.0 

DO  600  I-l.IM 
DO  600  J-l.JM 
DO  600  K“1.KM 
IF(E( I.J.X)-B)90.600.90 


90  L*I 


N»K 

NNN*0 

100  1H»XABSF(L1+L) 

IIC-XABSF(L2+M) 

IL*XABSF(L3+N) 

IF( IH)110.500.110 
110  IF{ IM-1H)500.120.120 
120  IF(  10130.500.130 
130  IF( JM-IK)500.140.140 
140  IF( IL)150.500.150 
150  IF<KM-IL)500.160.160 
160  1F<E( IH. IK. ILI-B 1 165.500. 165 
165  KH»XABSF(M1+L) 

KK«XABSF(M2+M) 

KL«XABSF<M3+N) 


1 


i 

f 
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167  IF(KH)170*500»170 

170  IF(  IM-ICH»500#18,0tl80 

180  IF(KK)190«500(190 

190  IF(JM-ICK I  500.200*200 

200  IF(KL)21Ui5U0.210 

210  IF(KM-KL)500.220.220 

220  IF(EIKH.KKiKL)-B>230.500>230 

230  RN‘RN+1.0 

TP»TP+E(  I  .J.(C)*E(  IH.U.IL)*E(ICH.KK.KL) 

500  NNN«NNN-fl 

GO  T0(510i520.510.530.510.540.510.600) .NNN 


510 

L»-L 

M«-M 

N«-N 

GO 

TO 

100 

520 

L«-L 

GO 

TO 

100 

530 

M»-M 

GO 

TO 

100 

540 

N»-N 

GO 

TO 

100 

600 

CONTINUE 

TP»TP*SF/RN+COR 

NR-RN 

TPC«TP/EEE 

I F ( XMODF ( JJJ  » 15 1 -1 1 620 .610 .620 
610  WRITE  OUTPUT  TAPE  JTP.llOO 

620  IF(MN)625.62S.621 

621  WRITE  OUTPUT  TAPE  JTP • 1300.M1.M2 .M3.N1 .Na.NS.NR »TP. TPC.COR 
PRINT  1300.M1.M2.M3.N1.N2.N3.NR.TP.TPC.COR 

GO  TO  9 

625  IFISENSE  SWITCH  1)626.627 

626  PRINT  1200<Ml*M2*M3iNl.N2.N3*NR.TP*TPC.COR 

627  Ll—Ll 
L2—L2 
L3»-L3 

WRITE  OUTPUT  TAPE  JTP . 1400.L1.L2 .L3 

WRITE  OUTPUT  TAPE  JTP * 1200.M1.M2 .H3.N1 *N2 .N3.NR *TP* TPC.COR 
GO  TO  9 
630  CALL  EXIT 

1100  F0RMAT{1H1.10X.75H  H  X  L  H  X  L  NO.CONTR. 

ID.  COSdNV.)  COR. TERM///) 

1200  FORMAT (IHO.IOX. 13.214. I 5.2I4.3X* 18. 2X *3( 3X.F10.5 ) ) 

1300  FORMAT(lH0*10X.I3.2I4*I5.2I4*3X*I8.2X.3(3X.E10.4) ) 

1400  FORMAT(1H0.10X*I3*2I4) 

END 


Tr4>le  Product  Summation,  Noncentrosymmetrlc 

CORTH  Rpe5§IOT?§8I5oI55T*  chwt. 

JTP-6 

_ KJPyj _ 

REWIND  KTP 

B«1«0E-21 _ 

READ  INPUT  TAPE  ITP.IOOO.NA.IM.JM.XM 

I  iMgjFORNAlLtULOJ _ 

i  DO  7S  I-l.IM 

I  DO  75  J»1.JM _ 

I  DO  79  R-l.KM 

_ 79  E(l.J.m«B _ 

r  READ  INPUT  TAPE  KTP.2000.N0IB 


i 


TRIP#PRO 


I 
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2000  FORMATdT)  _ _ _ 

00  8  NO-ltNOlB 

READ  INPUT  TAPE  KTP»2100»  1 1 J«IC»A _ 

2100  FORMAT(3I4tF20.9l 
IF(I)1»8»2 

1  1»-I 

2  1F(J)3»8»4  _  _ _ _ 

3  J—J 

4  IF(K»5f8t6  _  ,  _ 

9  R*~K 

6  E<I»J»K)»A  _ _ _ 

B  CONTINUE 

REWIND  KTP  _ _ 

AN>NA 

SF«SORTF(AN*»3)_/2.0 _ 

JJJ-0 

MN«0  _ 

C  START  NEW' TP' CALC. 

9  READ  INPUT  TAPE  I  TP. lOOOtMl tM2»M3*NU»d»N3 
JJJ-JJJ41  . . . 

Il»XABSF(Ml»  _ 

I2«XABSF(M2| 

I3«XABSF(M3)  _ _ _ _ 

J1>XABSF(N1) 

J2>XABSF(N2I 

J3*XABSP(N3) 

Ll»Ml+Nl  _ 

L2-M24N2' 

L3>M34N3 

Xl»XABSF<LH 

K2>XABSF(L2I 

X3«XABSF(L3» 

IF(Il»I24l3-»Jl»J2*J3H0tl0»90 _ 

10  MN-MN41 

_ EEE«lsP  _  _ 

CbR"0*b 

00  TO(99»690)»MW _ 

9b  iF(Ml)ll*20tll 

11  lF(N2)12t20tl2 _ 

12  lP(M3)13t20.13 

13  lFLAl>14«2a*-lA - 

14  lF(N2)19t20fl9 

19  IF(N3I17*20»17 - - - - - 

17  1F(L1)1B*20«18 

18  lF41.an*»30»lB - 

19  rF(L3)29»20t29 

20  CCC«1*0 - - 

COR-0.0 

00  TO  9»  - 

29  EEE>SORTF( (E( I 1*1 2*131+1. O)*(E(Jl»J2»J3)4l.0)*IE(Rl*K2tK3)4l.0t) 


CORa(ElIl.I3.I4l«F(Jl.J7.J4l«FIIC^  .K?  f -t-l  .(>  >  ^SORTF  f  f  N  I 

C  CALCULATE  TRIPLE  AVERAOES 

99  RN>0.a_ 

TP*0«0 

DO  600  iBl^ON. 

DO  600  J-1»JM 
_  DO  600  Fal.XM 

IF{En»J*K)*6)90*600t90 

90  L.-1  - . 

M>J 

NalL _ 

100 .  IH-XABSFILl+LJ _ 

IK*XABSP(L2+M} 

1L-XABSFLL3+I11. _ 

IF(tH)110t900»110 
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110  IF(IH-IH)S00fl20tl2O 

120  1F( IK)13O«S0O«130 

130  IF( JM-IK)5OO(140il4O 

140  IFdDlSOtSOOtlSO 

150  IF(i(H>IL)500tl60<160 

160  IF(E(IH*IK(IL)-Bll65t500tl65 

165  KH-XABSF(M1'«'L) 

KK>XABSF(M2+M| 

KL>XABSF(M3-t-NI 
167  IF(KH)17O*S0Oil70 
170  IF(IM-KH)500tie0iie0 
leo  lF(KK)190t500il90 
190  1F( JM-KK>500t200i200 
200  IF(KLI210t500t210 
210  IF(KM-XL)500>220i220 
220  IF(E<KHiKKiKL)-B)230i500*230 
230  RN-RN-t'l.O 

TP-TP'fEC  1iJ>K)*E(IH«IK>ILI»E(KH«XK«KL) 

500  NNN-NNN-fl 

GO  TO(510tS20«510«S30»510»540»510t600) tNNN 
510  L«-L 
M»-M 
N«-N 

GO  TO  100 
520  L»-L 

GO  TO  100 
530  M—M 

GO  TO  100 
540  M»-M 
N«-N 

GO  TO  100 
600  CONTINUE 

TP«TP*SF/RN+COR 

NR>RN 

TPC-TP/EEE 

IF(XMO0F(JJJt20) -1)620 #610*620 
610  WRITE  OUTPUT  TAPE  JTP.llOO 

620  IF(MN)625t625#621 

621  WRITE  OUTPUT  TAPE  JTP* 1300»M1.M2 .M3.N1 #N2 #N3.NR»TP.TPC#C0R 
PRINT  1300(Ml#M2#M3#Nl#N2#N3»NR#TP#TPCtC0R 

GO  TO  9 

625  WRITE  OUTPUT  TAPE  JTP# 1200#M1 tM2*M3>Nl >N2>N3*NR »TP# TPCtCOR 
PRINT  1200#M1#M2#H3#N1#N2#N3#NR#TP#TPC#COR 
GO  TO  9 
630  CALL  EXIT 

1100  FORMAT! 1H1#10X#75H  H  K  L  H  X  L  NO#CONTR#  TRIP. PRO 

10.  COSIINV.)  COR. TERM///) 

1200  FORMAT(lH0#10X#I3#2I4#I5#2I4#3X#ia#2X#3(3X#F10.5)  ) 

1300  FORMAT(1HO#10X#I3.2I4#I5#2I4#3X#I8#2X#3(3X#E10.4) ) 

END 


Structure  Factor  Calculation 

CSTRFAC  STRUCTURE  FACTOR  EXECUTIVE  PROGRAM 

DIMENSION  NAME(12) #X(3#500)#H(3#126)#  FI ( 20) #FA1 ( 20 ) #FB( 20) .FBI ( 

1  20)#FC(20)#  AA(3#3)#IA(500)#TB(500)#AF(500)#AF1(500)#BF(500 

2  )#6F1(500)#CF(500)#SI  (128).#E(128)#F(128)#A(128)#FA(128)#P(500)# 

3  GYP(500)#FX(128)#C(128)#SF(128) #EN 1 128 ) #B(  128 ) #0(128)#DF( 128)#DE( 

4  128)#CO(128)#DOI128)#DC(128)#DD(12B)#AO(128)#BO(128)#DA(128) tOBd 

5  28)#NMEd2)  #FJd28)#Aa(3#3)#SCF2d00) 
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50  HT*10 
IT*5 
JT»6 
LT  =  U 
REWIND  KT 
REWIND  LT 
NASP»0 
I»0 

READ  INPUT  TAPE  I T . 1000 . ( NAME ( 1)  .  I *l . 12 ) 

READ  INPUT  TAPE  I T t 1 100 • IC» I E«NR (NF t IPtLl B tLl NE 
READ  INPUT  TAPE  I T t 1200 t A1 • A2 tA3 <AL tBEtOA 
READ  INPUT  TAPE  I T t 1300 ,BT *XT .SCF » lOUT 
60  1«I+1 

READ  INPUT  TAPE  IT«1600,SCF2 ( I ) 

IF(SCF2( I  1)60.100.60 
100  NASP*NASP-*-l 

READ  INPUT  TAPE  I T . 1400 .FI ( NASP ) .FAI ( NASP 1 .FB ( NASP ) .FBI ( NASP 1 . 

1  FC(NASP) 

IFIFIINASP) 1 100.110.100 
110  K=0 
120  IC*K+1 

READ  INPUT  TAPE  I T  .  1500  .X  (  1  .K)  .X  ( 2  .K  1  .X(  3  .K )  .TB  ( 1C)  »  I A  (  K  ) 

IF(ABSF(X( l.K) )>ABSF(X(2.X) ll-ABSFIXI 3.K) ) >130.130.120 
130  X-X-1 
NC*X 

140  READ  INPUT  TAPE  I T . 1600 .T 1 .T2. T3 .U1 .U2.U3 

IF(ABSF(T1)+A9SF(T2)+ABSF(T3)+A9SF(U1)+ABSF(U2 )+ABSF(U3) >170.170. 

1  150 

150  DO  160  I«1.NC 
IC»X+1 

X( 1.K)*T1+U1*X( l.I ) 

X(2.X1=T2+U2*X(2.I) 

X(3.K)=T3+U3*X(3.I) 

IA(K)  =  IA( I ) 

160  TB(K)»TB(I) 

GO  TO  140 
170  SDF-O.O 
SF0*0.0 
SDF0=0.0 
DO  180  1=1. X 
J>IA(  I ) 

AF(n-Fl(J) 

AFK  I  )«FA1(  J) 

BF(I)=FB(J) 

BFK  I  )»FB1(  J) 

180  CF(I)=FC(J) 

IF(1E)190.195.190 
190  NN«1 

GO  TO  200 
195  NN*2 

200  CALL  RECIPI A1.A2.A3.AL.be. GA.ABI 1.1 ).AB( 2.2) .AB( 3.3) .AB( 2.3) .AB(1. 
1  3).AB(1.2)) 

DO  230  1=1.3 
DO  230  J*1.3 
IFt I-J)220f210.220 
210  AA(I>I)=  AB(I.I)»»2 
GO  TO  230 

220  AA(I.J)«  2.0*AB( I .J)*AB( I .1 )*AB( J.J) 

230  CONTINUE 
IF(NF)240.240.231 

231  DO  232  I>ltNF 

232  CALL  FLSXPD(XT) 

240  IF(LIB)242.243.242 

242  READ  INPUT  TAPE  KT.IOOO. ( NME< 1) . I>1 .12 ) 

243  GO  T0( 244.2401 ) .NN 
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244  IF(LINE)2403«2402<2403 

2402  READ  INPUT  TAPE  KT*2100«NR 
GO  TO  2401 

2403  READ  INPUT  TAPE  IT. 2100, NR 
2401  NO-0 

Nl-NR/128 

245  IF(N1-N0)470.250.260 
250  N2«NR-N1*128 

IF(N2»470»470,270 
260  N2-12e 

270  GO  70(280,310  ),NN 

280  IF(LINE)281,2e2,281 

281  READ  INPUT  TAPE  IT ,1700. ( H( 1 ,1 1 ,H(2 , 1 ) ,H( 3 , 1 ) ,E I  I ) , I >1 ,N2 I 
GO  TO  283 

282  READ  INPUT  TAPE  KT ,1700, ( H( 1,1) ,H( 2, I ) ,H( 3 , 1 ) ,E ( 1 1 , I *1 ,N2 ) 

283  DO  300  1-1, N2 

Ed)-  SORTF(E{I)+1.0)*SCF 
SI(  n«o.o 
00  290L-1,3 
00  290M-L,3 

290  SI(  n-H(L,I)*H(M,I)*AA(L,M)-«-SI(I) 

300  Sm)-SI(I)/4.0 
GO  TO  330 

310  IF(L1NE)311,312,311 

311  DO  2311  I-1,N2 

READ  INPUT  TAPE  IT, 1800,  H( 1,1 ) ,H( 2 , I ) ,H ( 3 , I ) ,F ( I ) , SCF.Sl ( I ) 
IF(SCF)2314, 2314, 2313 

2313  JK-INTG(SCF) 

2314  F(I)-S0RTF(F(n)»SCF2(JK) 

2311  CONTINUE 

GO  TO  313 

312  00  3311  I«1,N2 

READ  INPUT  TAPE  KT,1800,  H( 1 ,I ) ,H( 2 , I ) ,H ( 3 , I ) ,F ( I ) , SCF.Sl ( 1 ) 
IF(SCF)3314, 3314, 3313 

3313  JK-INTG(SCF) 

3314  Fd)-S0RTF(F(I))*SCF2(JK) 

3311  CONTINUE 

313  IF(SI(1) )314,314,317 

314  DO  316  I-1,N2 
SKD-O.O 

00  315  L-1,3 
DO  315  N-L,3 

315  SI(I)-H(L,1)*H(M#1)*AA(L,M)+SI(I) 

316  SI( I)-SI (I)/4.0 
GO  TO  330 

317  DO  320  I-1,N2 
320  SKI)-  S1(I)**2 
330  DO  450  I-1,N2 

A(I)-0.0 
SIGK-0.0 
FK(  n-o.o 

DO  350  J«1,K 

FA(  J)-AF(J)*EXPF(-AF1( J)*SI(I))+BF(J)*EXPF(-BFl(J)*Sni ) )*CF( J) 
SIGK-  SIGK'»'FA(J)*«2 

FX(  l)-FKd)+(FA(  J)*EXPF(-TB(J)*SIi(n)  )»*2 

P(J)-0.0 

DO  340  L-1,3 

340  P(J)-P(J)+H(L,I)*X(L,J) 

P(J)-P(J)*6.2831853 
GYP(J)-FA(J)*EXPF(-TB(  J)»SKI)) 

350  A(I)-A(I)4GYP(J)»C0SF(P(J)) 

FBTX-EXPF((BT*(SI(I))**XT)/2*0)/SQRTF(SIGK) 

FFAX>1,0/S0RTF(FK(I)) 

IF(KKT-1)3S1,352,353 
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351  FK(n»FFAX 
FJ(I)«FBTX 
GO  TO  355 

352  FK(n«FFAX 
FJ(  n»FFAX 
GO  TO  355 

353  FK(I)=FBTX 
FJ( 1 l-FBTX 

355  C(  n»A(  n*FK(  I  I 
1F( ICI370.360.370 
360  SF( I )-  A( I ) 

EN(n«  C(I) 

GO  TO  390 
370  B(n«0,0 

00  380  J-l.K 

380  Bdl-Bd)*  GYP( JI*SINF(P(Jn 
Dd)»Bdl*FK(n 
SF  d ) »S0RTF ( A ( d  **2+B { 1) **2 1 
EN(n«SFd)*FXd  ) 

390  GO  TO(400.410).NN 
400  F(  !)■  EdI/FJdl 
GO  TO  420 

410  Ed)»  Fd)*FJd| 

420  DFd)«SIGNF<Fd  I  .SF(  lll-SFd  ) 

SOF«SDF+ABSF(OFd  )  J 
SFO»SFO+F(n 

OEd  l=SlGNF(Ed  )  .ENd  d-EN(  n 
IF(F(  1  1  )422.423.422 

422  SDF0»S0F0+ABSF(DF(n ) 

423  IFdC)425.424.425 

424  F(  n»SIGNF(Fd  ).SF(I  )  ) 

E(n=SIGNF(E(n.EN(d) 

GO  TO  450 

425  IFdP)440.430.440 
430  TAR«Ed)/ABSF(EN(ni 

COd)»Cd)*TAR 
DOd)>Dd)*TAR 
DC(  n*SIGNF<C0(  I)  .C(  n  1-C(  I  I 
DOd  )  ■SIGNF(00(  I  I  .0(1 )  )-0(  n 
GO  TO  450 

440  RAT»Fd)/ABSF(SF(  n  ) 

AOd)»A(n*RAT 

BOd)>Bd)«RAT 

OAd)-SIGNF(AOd  ).A(n  )>A(n 
0Bd)>SIGNF(B0d  I.Bd))-B(  n 
450  CONTINUE 
460  IFdOUT)62t465.62 
62  IFdC)66. 65*66 

65  WRITE  OUTPUT  TAPE  LT .  1200.  (Hd .1  )  .H( 2 *  1)  .H( 3. 1  )  .F d  )  .Ed  )  .OF d  )  . 
10Ed)*I«l.N2) 

GO  TO  465 

66  IFdP)68.67*68 

67  WRITE  OUTPUT  TAPE  LT. 1200. ( Hd . I ) .H( 2 . 1 1 »H( 3. 1 > .CO (I ) .00(1 ) .DC (I ) * 
10D(I).I«1.N2) 

GO  TO  465 

68  WRITE  OUTPUT  TAPE  LT* 1200. ( H( 1. 1 ) .H( 2  *  I) *H( 3. 1 ) .AO ( I ) .BO( I ) .OA ( I ) * 
10Bd).I«UN2) 

465  CALL  OUTSFC(N2*IC.H.F*SF.A*B.E.EN.C*D*OF.NO.AO*BO.DA*OB*CO.OO.OC. 

1  DD.NAME.IP.OE) 

NO-NO-fl 
GO  TO  245 
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470  REWIND  KT 
END  FILE  LT 
REWIND  LT 
R-SDF/SFO 
Rl-SDFO/SFO 

WRITE  OUTPUT  TAPE  JT . 1900.SFO*SDF,SDFO»RfRl 
READ  INPUT  TAPE  ITt2100«njJ 
IF( IIJJ)50<480i50 
480  PRINT  2000 
CALL  EXIT 
1000  FORMAT! 12A6) 

1100  F0RMAT(7I101 
1200  FORMAT(7F10.4) 

1300  FORMAT(3E10.4»2I10I 
1400  F0RMAT(5F8.4) 

1500  FORMAT(4E1Q.4.I10) 

1600  F0RMAT!6E10,41 
1700  FORMAT(3F4.0.F20.5) 

1800  FORMAT(4F9.2 .9XtF9«2»F9.6» 

1900  FORMAT! 1H1.9X.12HSIGMA  FIOl  »*F15.3»eX . 15HS1GMA  /DEL  F/  ■tF10.3.8X 
1»15HS1GMA  /del  F/  =.F10.3»1H*///24X,2!3HR  « »F9. 5 »9X ) . 1H*////20X t64 
2H*  THESE  QUANTITIES  00  NOT  INCLUDE  CONTRIBUTIONS  FROM  REFLECTIONS/ 
320X.33H  OBSERVED  WITH  ZERO  INTENSITIES.) 

2000  F0RMAT!1H1.10X«21HTHIS  JOB  IS  FINISHED.///////////) 

2100  FORMATII7) 

END 


Subroutine  OUTSFC 

SUBROUTINE  OUTSFC I N2 . 1 C.H.F.SF. A.B»E»EN.C »D.DF »N0» AO.BO»DA,DB»CO» 

1  DO.DCiDD.NAME.IP.OE) 

DIMENSION  H!3il28)«Fll2e).A!12B) .B!128)>E!128)»EN!128)»C!12B)>D!12 

1  8).  DF!128)»AO!128)»BO!128)»DA!128)»DBI128)»COI128)»00!128) 

2  .001128) *00! 128)  .NAME  1 12 ) .DE ! 128 ) .M! 3. 128 ) .SF ! 128  ) 

JT=6 

N4»6*N0 
DO  50  I«1.N2 

DO  50  J»l,3 

50  M!J.I)>INTG!HIJiI  )) 

IF!128-N2)110.110.100 
100  N1-N2/22 

N3»N2-N1*22 
IF!Nl)165«165il20 
110  Nl«5 
N3-18 

120  DO  160  I-1«N1 
N4«N4+1 
J1«!I-1)*22+1 
J2-J1+21 

WRITE  OUTPUT  TAPE  JT . 1000 ,! NAME! X) »X» 1 » 12 )  .N4 
IF!  10125.150.125 
125  IF! IP)140.130.140 
130  WRITE  OUTPUT  TAPE  JT.1200 

WRITE  OUTPUT  TAPE  JT . 1400. ! M ! 1 » J ) .M! 2 » J ) .M! 3. J ) »F ! J ) ,SF ! J) .E! J ) . 

1  EN!J).CO!J)*DO! J).C! J)*0!J).DC! J).DD!J).DE! J).J>J1.J2) 

GO  TO  160 

140  WRITE  OUTPUT  TAPE  JT.llOO 

WRITE  OUTPUT  TAPE  JT . 1300. ! M! 1 » J ) .M! 2 1 J ) .M! 3. J ) »E ! J )  .EN ! J ) .F ! J ) . 

1  SF! J).AO!J) .BO! J).A! J)*B!J)»DA! J)«DB! J) .DF! J).J-J1*J2) 

60  TO  160 

150  WRITE  OUTPUT  TAPE  JT.1500 

WRITE  OUTPUT  TAPE  JT. 1600. IM! 1 » J ) .M! 2 . J ) »M! 3. J ) .F ! J) .SF ! J) .E ! J ) . 

1  EN!J).DF!J)*DE!J).J«J1.J2) 

160  CONTINUE 
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165  IF(N3)210»210.170 
170  J1«N1*22+1 

N4»N4+'l 

WRITE  OUTPUT  TAPE  JT t 1000 » ( NAME( K I *K«1» 12 ) f N4 

IF(1C)175.2C0.175 
175  IF(lP)190.180il90 
180  WRITE  OUTPUT  TAPE  JTtl200 

WRITE  OUTPUT  TAPE  JT . 1400 1 ( M ( 1  . J ) »MI 2 ♦ J) .M ( 3 * J ) *F ( J ) ,SF ( J ) tE ( J ) * 
1  EN(J) tCO( J)«DO( J)«Ci JltD( JltOC(J)«DO(J) •DE(JI*J>JliN2l 

GO  TO  210 

190  WRITE  OUTPUT  TAPE  JT.llOO 

WRITE  OUTPUT  TAPE  JT 1 1300 « (M ( 1 • J I (MI  2 tJ ) tMI 3 • J ) •£ ( J ) lEN ( J ) fF ( J I t 
1  SF( Jl tAOI J) tBOI J) tA( J) tBl J) tOA( Jl fOBIJ) «DF( J) «J*JltN2) 

GO  TO  210 

200  WRITE  OUTPUT  TAPE  JTtlSOO 

WRITE  OUTPUT  TAPE  JT . 1600 1 ( M ( 1 . J ) .MI 2 . J ) tM ( 3 1 J ) tF ( J ) ,SF ( J ) ,E ( J ) , 
1  EN(J)i0F(J)t0E( J)>J«JliN2) 

210  RETURN 

1000  FORMAT(1H1«12A6//100X»4HPAGE»I4/) 

1100  FORMAT(120H0  H  K 
1A(0)  8(0) 

2) 

1200  F0RMAT(120H0  H  K 
1C(0)  0(0) 

2) 

1300  FORMAT(lHOt313f2F10.4*9F10.3) 

1400  FORMAT) 1H0» 313 •2F10.3«9F10. 4) 

1500  FORMAT(86HO  H  K  L  F(0) 

1  E(C)  OEU  F  DEL  E  /) 

1600  F0RMAT(2H0  f3I4*2F12.3*2F12*4»F12.3»F12.4) 

END 


L 

A(C) 

E(0) 

B(C) 

E(C) 

DEL 

A 

F(0) 

DEL 

B 

F(C) 

DEL 

L 

C(C1 

F(0) 

D(C) 

F(C) 

DEL 

C 

E(0) 

DEL 

D 

EIC) 

DEL 

F(C) 


E(0) 


Interatomic  Distance  and  Angle 

CBNDANG  BONO  DISTANCE  AND  ANGLE  PROGRAM 

DIMENSION  NAME( 12) tN(500) tX(500).Y(500) . 

1Z(500),NS1(200) tNS2(200) tSV ( 200 ) *N1 ( 2900 ) f N2 ( 2900 ) * 

2D( 2900 ),E1( 2900 ) tE2( 2900) .E3( 2900) ,ANG(2000) * 

3M1(20C0) »M2( 2000)  .M3 (2000) 

COMMON  NltN2>DtEl.E2tE3 

EQUI VALENCE ( Ml. N). (M2. X) . ( M3 . Y ) . ( ANG.Z) 

COSOF(X)=COSF(3.14159*X/180.) 

1TP=5 

JTP*6 

XTP»9 

REWIND  XTP 

READ  INPUT  TAPE  I  TP. 1000 . ( NAME ( 1 1 ) » 1 1«1 » 12  ) 

READ  INPUT  TAPE  1  TP. 1200 • A.B.C .AL.BE tGA 

READ  INPUT  TAPE  I  TP . 1500 .BMX . AMX 

READ  INPUT  TAPE  ITP.llOO.XINC.YINC.ZINC 

BMX»BMX**2 

CAL«COSDF(AL) 

CBE»COSDF(BE) 

CGA*CbS0F(GA) 

I»0 

IF { XI NC+YINC+Z INC  1107.100.107 
100  I«I+1 

READ  INPUT  TAPE  ITP.llOO .X( I ) . Y { I ) .Z( I ) 

105  N(n«I-l 

IF(ABSF(X(I  )  )-«-ABSF(Y(  I ) ) 4ABSF ( Z(  1 ) ) ) ll'O .  110. 100 
107  CALL  iNTERPd.N.X.Y.Z.ITP.XINC.YlNC.ZINC. A.B.C. CAL. CBE.CGA) 


02693 

02693 

02693 


02693 

02693 

02693 

02693 

02693 

02693 
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110  JJ-0  02693 

KK-0  02693 

I- I-l 

II- I  02693 

115  READ  INPUT  TAPE  ITPil200,Tl»T2»T3»Ul»U2«U3 

lF(ABSF(Tl)+ABSF(T2)+ABSF(T3)+ABSF<Ul)4.ABSF(U2)+ABSF(U3))140tl40*  02691 
1  120  02691 

120  DO  130  J>1»I  02693 

Il-Il+l  02693 

N(I1)>I1-1  02693 

X(  Il)«Tl-«-Ul*X(  J)  02691 

Y( I1)»T2+U2*Y< J1  02691 

130  Z(I1)>T3>U3*Z(J)  02691 

60  TO  115  02693 


140  DO  145  J-1«I1 

IF (XMODFIJ. 50 )-l)142» 141,142 

141  WRITE  OUTPUT  TAPE  JTP, 1800, ( NAME (K ) ,X-1 , 12 ) 

142  LlaN(J)/I 
L2-N( J)-L1*I+1 

145  WRITE  OUTPUT  TAPE  JTP » 1700,L1»L2 »X< J ) » Y ( J ) ,2 ( J ) 


150  D0190  J-1,I1  02691 

I2-J-«-l  02693 

DO  190  K-I2,I1  02693 

IF(N(J)/1-N(K) /I) 157,155,157 

155  IF(N( J)/ni90,157,190 

157  Vl'X(Jl-X(IO  02693 

V2*Y(J)-Y(X1  02693 

V3»2(  Jl-Z(IC)  02693 

VEC>ABSF(00TPR0( V1,V2,V3,V1,V2,V3,A,B,  02693 

lC,CAL,CBE,C6An  02693 

IFIVEC-  BMX)160,160,190 

160  IF(VEC-.49)170,180,180  02693 

170  JJ«JJ+1  02693 

NS1(JJ)-N(J)  02693 

NS2(JJ)«N(K)  02693 

SV( JJI-SQRTFIVFf'  02693 

GO  TO  190  02693 

180  XKoKK'fl  02693 

N1(KK)>N(J)  02693 

N2(KK)-N(X)  02693 

D(KK)«SQRTF(VEC)  02693 

E1(KK)-V1  02693 

E2(KK)«V2  02693 

E3(XX)-V3  02693 

I F(2900-XK) 330,330,190 

190  CONTINUE  02693 

LL>0  02693 

NM«0 

00  322  J«il,KK 
IF(D(J)*AMX) 195, 195,322 

195  XKl>J-t-l  02693 

DO  320  K-KK1,KK  02693 

IF(D(K)-AMX 1200,200, 320 

200  IF(N1(J)-N1(K)|220,210,220  02693 

210  M>1.  02693 

ZD-1,0 
60  TO  276 

220  IF(N1(J)-N2(K))240,230,240  02693 

230  M-2  02693 

ZD— 1.0 
GO  TO  276 

240  lF(N2(J)-Nl(ICn260,290,260  02693 

250  M-3  02693 

ZD— 1,0 
60  TO  276 
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260  IP(N2(J)-N2(Kt >320t270t320 
270  M-4 

ZD-1*0 

276  LL«LL*1 

C0ANG>ZD«  OOTPROtEK J)tE2(JI tE3IJ) 

1E2(K) *E3(K) •AtB*C«CALtCBE*CGA)/(0( J|*0(Kn 
ANGC0«1.0>ABSF(C0ANG) 

IF(ANGCO»2277t277f277 

2277  WRITE  OUTPUT  TAPE  JTPi  1600»N1 ( J) »N2 ( J ) tNl ( H ) »N2 ( K ) tANGCO 
COANG-1.0 

277  ANG(Lt»»90, 00-57. 295780*ARCSIN(COANG» 

278  GO  T0(280t290«300«310)  iM 
280  M1(LL)>N2(JI 

M2(LL)-N1( J) 

M3(LLI-N2(K) 

GO  TO  315 
290  H1(LL)-N2(J) 

M2(LLI«N1( J) 

H3(LL)>N1(KI 
GO  TO  315 
300  MllLLt-NUJ) 

M2(LL)>N2( J) 

H3(LL)«N2(K) 

GO  TO  315 
310  M1(LL)«N1(J) 

H2(LL)-N2( J) 

M3(LL)«N1(K) 

315  IF(2000-I.LI360t360t320 

360  WRITE  TAPE  KTP* ( Ml (MN ) •M2 (MN ) tM3 (MN ) »ANG(MN ) >MNsl»2000 ) 
LL«0 
NM«NM4l 
320  CONTINUE 

322  CONTINUE 
IF{LL)323. 323. 1323 

1323  WRITE  TAPE  KTP • ( Ml (MN ) •M2 ( MNl •M3 <MN) (ANGIMN) •MN>1.LL ) 

323  CALL  OUTBND(I.JJ*KK.LLtNSl.NS2.SV. 
lNliN2.N3.0.El.E2.E3.ANG«Ml«M2.M3. 

2NAME.JTP.NM>KTP) 

325  CALL  EXIT 

330  BMX«<S0RTF(BMX)-.25)**2 
IF(BMX-4.0 1335.335. 370 


02693 

02693 


02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 


02693 


02693 

02693 


335  IF(I1/I-1»340.340.350 
340  WRITE  OUTPUT  TAPE  JTP.1300 
GO  TO  325 

350  WRITE  OUTPUT  t/vpe  JTP.1400 
I1>I 
370  KK«0 
JJ»0 

GO  TO  150 
1000  FORMAT! 12A6) 

1100  F0RMAT(3E10.4) 

1200  F0RMAT(6E10.4I 
1300  FORMAT! 1H1.8X.71HM0RE  THAN  2900  2  A*  OR  LESS  BOND  DISTANCES  HAV 

IE  BEEN  FOUND  AMONG  THE///9X.71HUNTRANSLATED  SET  OF  ATOMS.  THIS  PR02691 
20BLEM  IS  TOO  LARGE.  NO  RESULTS  WILL///9X'.10HBE  ISSUED.)  02691 

1400  FORMATI1H1.8X.71HMORE  THAN  2900  2  A.  OR  LESS  BOND  DISTANCES  HAV 

IE  BEEN  FOUND.  VECTORS///9Xi58HWlLL  BE  CONSIDERED  FOR  THE  UNTRANSL02691 
2ATED  SET  OF  ATOMS  ONLY.)  02691 

1500  F0RMATI2E10.4) 

1600  FORMAT! 10X.4I4.E12.5) 

1700  F0RMAT!11X.I2.1H-.I2.3F10.5) 

1800  FORMAT!1H1.10X.12A6///12X.4HATOM.6X.1HX.9X.1HY.9X.1HZ/) 


02691 


02691 

02693 

02693 

02693 


END 
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Subroutine  INTERP 


SUBROUTINE  I NTERP (I (N tX tY tZ • ITP.XI NC t YI NC tZ I NC • A«B tC tCALtCBEtCGA) 
DIMENSION  X(500)  fYlSOOl  tZISOOItNCSOO)  iP(7)  fEOf  3)>G(3«4)  tHO) 

E( 

E( 1«2)»B*CGA*YINC/(X1NC*AJ 
E(1«3)=C*CBE*ZINC/(X1NC»A» 

E ( 2  *  1 1 =A*C6A*X I NC/ ( Y 1 NC*B » 

E(2«2)=1.0 

E(2t3)*C*CAL*ZINC/(YINC*B) 

El 3»1 1«A*CBE*X1NC/(ZINC*C) 

E(3»2)=B*CAL*YINC/(ZINC*C) 

E(3«3)«1.0 
100  I«I+1 

READ  INPUT  TAPE  I  TP t lOOO.X ( H *Y( 1 ) *2 1 1) 

IFIABSFIXn ) )+ABSF(Y( 1 D+ABSFI Z ( I ) ) >1 30 . 130. 1 10 
110  READ  INPUT  TAPE  I  TP. 1000 . IP (J) . J= 1 . 7 ) 

L=0 

DO  120  J=1.3 
L  =  L+2 

G ( J .4 ) »L0GF I P ( 1 ) /P ( L I  ) /LOGF ( P (  1 ) /P I L+1 ) ) 

G(J.4)«(G(  J.4I-1.0)/(2«0*(G(  J.4)-«-l*0)  ) 

DO  120  K=1.3 
120  G(J.K)«E(J.K) 

CALL  MATSIG.H.3.11 
X(  n»X(  1  )+H(  1)*XINC 
YII)-YII)+H(2)»Y1NC 
21  I l*Z( I I+H(3)*2INC 
Nm»I-l 
GO  TO  100 
130  RETURN 
1000  FORMAT(7E10.4) 

END 


Subroutine  OUTBND 


SUBROUTINE  OUTBNOI I . J J.XK .LM.NSl . 

1NS2.SV.N1.N2.N3.0.E1.E2.E3.ANG.M1. 

2M2iM3.NAME.JTP.NM.XTP) 

C  OUTPUT  S.R.  FOR  BOND  DISTANCE  AND  ANGLE  PROGRAM  02693 

DIMENSION  NSl(20a).NS2(200)»SV(200»  .  02693 

1N1(3000).N2(3000) .013000) »E1( 3000) tEZ 1 3000). 

2E3I3000) .ANGI 2000) .Ml (2000) .MZI 2000) . 

3M3(2000).NAME(12I  02693 

REWIND  KIP 

IJK»1  02693 

50  Il>KK/44  02693 

IF(KK)85i85.75  02693 

75  IF(Il)120.12OtlOO  02693 

85  WRITE  OUTPUT  TAPE  JTP.8000 

GO  TO  250  02693 

100  DO  110  J-l.Il  02693 

WRITE  OUTPUT  TAPE  JTP . 1000. (NAME IK ) .K=l . 12 ) . J 

13*J*44-22  02693 

I2-I3-21 

DO  110  L*I2.I3  02693 

L1«N1(L)/1  02693 

M11«N1(L)-L1*I+1  02693 

L2>N2(L)/I  02693 

M22«N2IL)-L2*I+1  02693 

LLl»Nl(L+22)/I  026'93 

MMl»Nl(L+22)-LLl*I  +1  02693 

LL2«N2(L-^22)/I  02693 

MM2>N2(L'«-22)>LL2»I  *l  02693 
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110  WRITE  OUTPUT  TAPE  JTP»2000tL 1 tMl 1 .L2.M22 *D( L) tLll tMMl .LL2 *MM2 i 
lD(L+22) 

IP(ICK-Il*44)170tl70*120 

120  14-11*44+1 
I5-KK-I4-21 
IF(I5)121>122«122 

121  15-0 

122  I6-I4+I5-1 

17- I6+1 

18- KIC-I5 
19*11+1 

WRITE  OUTPUT  TAPE  JTP. lOOOt ( NAME (K » tK-l ♦ 12 ) *  19 
IF(I5)150»150tl30 
130  00140  L-I4fl6 
L1-N1(L)/1 
M11-N1(L»-L1*I  +1 
L2-N2a)/I 
M22»N2(L)-L2*I+1 
LLl*Nl(L+22)/I 
MMl*Nl(L+22)-LLl*I  +1 
LL2*N2(L+22I/I 
MM2-N2(L+22)-LL2*I  +1 

140  WRITE  OUTPUT  TAPE  JTP i2000.Ll .Mil tL2 .M22 *0 ( L ) »LLl *MM1 . 

lLL2(MM2tO(L+22  ) 

150  DO  160  L*I7.I8 
L1-N1(L)/I 
M11*N1<L)-L1*I  +1 
1.2»N2(L)/I 
M22*N2(L)-L2*I  +1 

160  WRITE  OUTPUT  TAPE  JTP * 3000*L1 tMll »L2 »M22 *D ( L) 

170  GO  TO(175»270)*IJIC 

175  IF(NM)250tll76tll75 

1175  LL»2000 

GO  TO  1177 

1176  LL-LM 
IF(LLI250t250tll77 

1177  NM*NM-1 

READ  TAPE  KTP t ( Ml ( MN I *M2 ( MN I iM3( MN ) tANGI MN ) »MN- 1 »LL) 

Jl-LL/44 

176  IF(Jl)200t200»180 
180  DO  190  K*1,J1 

WRITE  OUTPUT  TAPE  JTP *4000. (NAME ( L) »L»1 » 12 ) 1 1C 

J3=IC*44-22 

J2-J3-21 

DO  190  N-J2.J3 

U1»M1(N)/I 

L2-M1 (N)-L1*I  +1 

L3*M2(N)/I 

L4*M2 (N1-L3»I  +1 

L5»M3(N)/I 

L6*M3(N)-L5*I  +1 

MLl*Ml(N+22)/I 

ML2»Ml(N+22)-MLl*I  +1 

ML3*M2(N+22)/I 

ML4»M2(N+22)-ML3*I  +1 

ML5*M3(N+22)/I 

ML6-M3(N+22)-ML5*I  +1 

190  WRITE  OUTPUT  TAPE  JTP. 5000.L1 »L2 .L3 .L4»L5*L6 »ANG( N ) * 
lMLl.ML2.ML3.ML4.ML5.ML6«ANG(N+22) 

IF(LL-Jl*44)175il75.200 

200  J4»J1*44+1 
J5-LL-J4-21 
•IF(J5)20l.20?,202 

201  J5-0 
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02693 

02693 

02693 

02693 

02693 

02693 


02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 


02693 


02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 

02693 
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202  J6«J4+J5-1 
J7-J6+1 
J8«LL-J5 

J9«J1+1  02693 

WRITE  OUTPUT  TAPE  JTP .4000 • (NAME (L ) *L«1 *  12 ) .J9 

IF(J3)230.230.210  02693 

210  DO  220  N-J4.J6  02693 

L1*M1(N)/1  02693 

L2«M1(N)-L1*I  +1  02693 

L3«M2(N)/1  02693 

L4-M2(N)-L3*I  +1  02693 

L5>M3(N)/I  02693 

L6«M3(NI-L5*I  +1  02693 

MLl-Ml(N-»'22)/!  02693 

ML2«Ml(N-*-221-MLl*I  -t-l  02693 

ML3-M2(N+221/I  02693 

ML4>M2(N-*-22)-ML3*I  *1  02693 

ML5>M3(N-*-22)/I  02693 

ML6«M3(N+22)-ML5*l  +1  02693 

220  WRITE  OUTPUT  TAPE  JTP • 5000.L1.L2 >L3 •L4»L3.L6.ANG( N )  • 

lMLl.ML2.ML3.ML4.ML5.ML6.ANG(N-f22)  02693 

230  DO  240  N-J7.J8  02693 

L1=M1(N)/I  02693 

L2»M1(N)-L1*I  +1  02693 

L3>M2(N)/I  02693 

L4«M2tN)-L3*I  +1  02693 

U5»M3(N1/I  02693 

L6»M3(N)-L5*I  *l  02693 

240  WRITE  OUTPUT  TAPE  JTP •6000.L1.L2 .L3.L4.t$.L6.ANG( N ) 

GO  TO  175 

250  IF(JJ)265<26S»2SS  02693 

255  KK«JJ  02693 

DO  260  J>1«JJ  02693 

N1(J)>NS1(J)  02693 

N2(J)>NS2(J)  02693 

260  D(J)«SV(J)  02693 

WRITE  OUTPUT  TAPE  JTP. 7000 

IJK-2  02693 

GO  TO  50 

265  WRITE  OUTPUT  TAPE  JTP. 9000 
270  REWIND  KTP 
RETURN 


lOOO  FORMAT! 1H1.6X.30NBONO  DISTANCES! ANGSTROMS)  FOR  .12A6.7H  PAGE  .12/02693 
1//19X.24HAT0M  ATOM  DISTANCE. 27X .24HAT0M  ATOM  DI STANCE ) 02693 

2000  FORMAT!1H0.17X»I2.1H-.I2.3X.I2.1H-.I2.4X.F8.5.26X.I2.1H>.I2.3X.I2> 
11H-.I2.4X.F8.5I 

3000  FORMATIIHO,17X»I2.1H-.I2.3X.I2.1H-.I2.4X.F8.5) 

4000  FORMAT! 1H1.6X.30HBOND  ANGLES  (DEGREES)  FOR  .12A6.7H  PAGE  .12/02693 

1//2!11X.32HAT0M  ATOM  ATOM  ANGLE. IIX))  02693 

5000  FORMAT!1H0.2!9X.I2.1H-.I2.3X.I2.1H-.I2.3X.I2.1H-.I2.4X.F9.4.10X) )  02693 
6000  FORMAT!1H0.9X>I2.1H-.I2.3X.I2.1H-.I2.3X.I2.1H-.I2.4X.F9.4)  02693 

7000  FORMAT! 1H1.10X.53HTHE  FOLLOWING  BOND  DISTANCES  ARE  ALL  LESS  THAN  002693 
1.7  A.//11X.59H80ND  ANGLES  INVOLVING  ThESE  BONOS  HAVE  NOT  BEEN  CALC02693 
lULATED.)  02693 

8000  FORMAT! 1H1.10X.64HN0  BOND  DISTANCES  BETWEEN  0.7  AND  4.0  ANGSTROMS  02693 
IHAVE  BEEN  FOUND.)  02693 

9000  FORMAT! 1H1.10X.56HN0  BOND  DISTANCES  LESS  THAN  0.7  ANGSTROMS  HAVE  B02693 
lEEN  FOUND.)  02693 

end  02693 
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Least-Squares  Line  and  Plane  Fitter 


clsopln  least-squares  line  and  plane  fitter 

DIMENSION  ZZ(12)  tG(3i3)(2Y(  lOUXOtSOItWXSO)  *  WXXS  ( 3>3  )  iWT  (  50 )  » 
1ADJ( 3*3) t  A(3i3) tGl(3«3l •XHK(3) <B(3«3)>YMK(3) •VM(3»3) (VMI (3>3) i 
2VMJ(3f3)>  C(3t3)<  01 1  50 ) tOEL ( 50 1 •OEL2 ( 50 ) tXBAR ( 3 ) • VMX(3 *50 1 > 
3P(3)*UMK(3*50)  •BV(3) 

4*Z12(3) 

ITP*5 

JTP-6 

Z34*1.0E-12 

MM*10 

READ  INPUT  TAPE  I  TP . 1000 , ( ZZ ( IJ . 1*1 .12 1 
READ  nPUT  TAPE  I  TP. 1003. A1 . A2 .A3.AL.BE.GA 
WRITE  OUTPUT  TAPE  JTP . 1020 . ( ZZ ( I ) . 1*1 .12 ) 

CALL  RECIPI Al.A2.A3.AL.BE.GA.Bl.B2^.B3*COSA*COSB.COSC) 

C  SET  UP  MATRIX  G 
G(1*1)*B1««2 
G(1.2)»B1*B2*C0SC 
G( 1.3)=B1*B3*C0SB 
GI2*1)*G(1*21 
G(2*2)*B2**2 
G(2.3)»B2*B3*COSA 
G(3*1)*G(1.3) 

G(3.21*G(2.3) 

G(3»3)»B3**2 

G1(1.1)*A1»»2 

G1(1.21*A1*  A2*COSF(GA*3. 14159/180,0) 

GKl  *3 1*A1*A3*C0SF(BE»3. 14159/180,0) 

G1(2*1)«G1(1,2) 

G1(2,2)*A2**2 

Gl( 2.3 )*A2*A3*COSF<AL*3, 14159/180,0) 

G1(3.1)*G1( 1*3) 

G1(3*2)-G1(2.3) 

G1(3.3)*A3**2 

READ  INPUT  TAPE  I  TP. 1002 *N *NO.L. (ZY ( I ) . I *1 *10 ) 

10  DO  20  J»1.N 

20  READ  INPUT  TAPE  ITP.1003,  X ( 1 . J) ,X (2 . J) .X( 3. J ) .WT ( J ) 

C  COORDINATES  OF  CENTROID 
DO  25  I-l.N 
IFIWTII) )25.22.25 
22  WT(I)>1,0 
25  CONTINUE 
DO  30  1*1.3 
WXS(I)*0. 

DO  30  J*1.N 

30  WXS(  I)*WXSn)+WT(J)*X(  I.J) 

WS*0. 

DO  40  1*1. N 
40  WS*WS-l-WT(  I ) 

DO  50  I«l,3 
50  XBARI1)*WXS(I)/WS 
C  SET  UP  MATRIX  A 
DO  60  1*1,3 
DO  60  J*1.3 
WXXS(I*J)*0. 

DO  55  K*1.N 

55  WXXSI I.J)*WXXS( 1 *J)-fWTIK)*X( I.K)»X(J.K) 

60  A(I.J)«WXXS(I*J)-XBAR(I)*XBAR(J)*WS 
IF(L)65*65.61 

C  EVALUATE  MATRIX  B  FOR  LINE 

61  CALL  MTXMUL(3.3*3*G1*A*B) 

GO  TO  75 

C  SET  UP  ADJOINT  OF  A 


64 


NAVAL  RESEARCH  LABORATORY 


65  ADJ( l»l»*A(2t2)*A{3»3)-A<2*3»*A(3»2) 

ADJ(2.1)»AI  3.1.)*A(2.3)-A<2.1)*At3f3) 
ADJ(3tl)-A(2tl)*A(3t2)-A{3.1)*A(2»2) 

ADJ( l»2)-A(3.2)*A(1.31-A(lt2)*A(3.3) 

ADJ(2*2I»A(  ltl»*A(3.3)-A(3.1»*Aa»3) 
ADJ(3»2)»A(3tl)*A(l»2)-A(l.l)*A(3.2l 
ADJ( 1»3)»A( 1.2)*A(2.31-A(1.3)*A(2»2) 
ADJ(2.3)»A(2tl)*A(1.3)-A(l*l)*A(2»3) 

ADJ(3.3|»A( ltll*A(2«2)-A<2.1)*A(l*2» 

C  EVALUATE  DETERMINANT  OF  A 

CALL  MTXMUL( 3*3t3iAOJ»AtC> 

DETNA=0, 

DO  70  1-1.3 
70  DETNA-DETNA+Cd.n 
DETNA-DETNA/3.0 

C  EVALUATE  MATRIX  B  FOR  PLANE 

CALL  MTXMUL(3.3»3.ADJ»G.B) 

75  DO  80  1*1.3 
80  BV(I)=0. 

DO  90  1-1.3 
DO  90  J-1.3 

90  BV(n=BV(I)+B(J.I)**2 
BI6BV=BV(1) 

sr'uo 

IF(BV( n-BIGBV) 110.110.100 
100  KK-I 

BIGBV=BV( I » 

110  CONTINUE 

VMN  =  S(JRTF(BIGBV) 

DO  120  1-1.3 
120  VM(I.1I-B(I»KK)/VMN 
NNN-0 

C  EVALUATE  VM  BY  ITERATION 

125  CALL  MTXMUL(3.3.1«B.VM»VMI) 

NNN-NNN+1 

IF(MM-NNN)3.3.4 

3  MM-MM+10 
234-234*100. 

4  IF(40-NNN)1.1.2 

1  WRITE  OUTPUT  TAPE  JTP.1021.NO 
GO  TO  235 

2  IJ»1 

VMN-SORTF(VMI ( 1.1 ) **2+VM I ( 2. 1)**2+VMI (3.1)**2 I 
DO  126  1-1.3 

126  VMI(I.1)-VMI(I.1)/VMN 
DO  140'l-1.3 

212  (  n»VM(  I  .1)/VMI  (1.11-1,0 

IF(ABSF  ((VM(I,ll/VMI(I»in-1.0)-234  )140»140.130 

130  IJ=2 

GO  TO  150 
140  CONTINUE 
150  DO  160  1=1,3 

VMd.D-VMKI.l) 

160  VM(1.I)-VM( I.l) 

GO  T0d70,125),IJ 
C  NORMAL12E  VECTOR  VM 

170  CALL  MTXMULd.3.3.VM.G.VMI  ( 

CALL  MTXMUL(1,3.1,VMI.VM,VMJ) 

ORM-SORTF(VMJ(l,l) I 
DO  180  1*1,3 
180  VMd.l)-VMd.l)/ORM 
IF(L)185.185.500 
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500  DO  510  1*1,3 
P(  11*0.0 
00  510  J*l,3 

510  PI  n*VM(J,l)*G(J,n+PI  1) 

DL2«0,0 
DO  530  1*1, N 
DEL2(n*0.0 
DO  520  J*l,3 
DO  520  K«l,3 

520  DEL2I I l*0EL2(I )> I XBAR IJ )-X (J , I ) ) •IXBARI  K.)-X ( K,  I  ))*(G1(J,X)- 
1VM(J,1)*VM(K,11 ) 

530  DL2*0L2+DEL2n  ) 

WRITE  OUTPUT  TAPE  JTP , 1025 ,NO,N, IZY ( I ) , I *1 , 10 ) ,XBAR ( 1 ) ,P 1 1 ) ,XBAR( 2 
II  ,P(2),XBAR(3),P(3I,(  I,WT(  n,DEL2(n*I>l*N) 

SATN*N-2 

IF(SATN)55{),550,540 
540  STN0*SQRTF(0L2/SATN) 

WRITE  OUTPUT  TAPE  JTP, 1032 ,STN0 
550  DO  560  1*1,3 
560  UMK(I,NO)*VMII,ll 
GO  TO  235 

C  EVALUATE  EIGEN  VALUE  FOR  PLANE 
185  CALL  MTXMUL(3,3,1,B,VM,VMJ| 

0NEG>0. 

DO  190  1*1,3 

190  0NEG=DNEG+VMJ(I,1)/VM(1,1) 

EIGEN=3,0*0ETNA/DNEG 

C  EVALUATE  DIPLANE  TO  ORIGIN  DISTANCE) 

0»0, 

DO  200  1*1,3 
200  0*0+VMI 1 ,1 )*X8AR ( 1 1 

C  EVALUATE  D, DELTA  D,IDELTA  01**2  FOR  EACH  POINT 
DL2*0. 

DO  220  1*1, N 
DI  ( 1 1=0, 

DO  210  J*l,3 

210  Oim=OIII)+VM(J,ll*XIJ,I) 

DELI  I  1=01 1  I  1-0 
DEL2I  n*DEL(  1 1**2 
220  DL2*DL2+WTIII*DEL2II) 

DO  230  1*1,3 

230  VMKII,NOI»VMII,l) 

C  OUTPUT  ONE 

WRITE  OUTPUT  TAPE  JTP , 1005 ,NO,N, IZY 1 1) , I »l, 10 1 , VM 1 1 , 1 1 , VMI 2 , 1 ) ,VMI 
13,1),D,XBARI1|,XBARI2I,XBARI3),I I ,WT I  1 1 ,DI 1 1 1 ,OEL 1 1 1 ,DEL2 1 I),I*1,N 
2) 

WRITE  OUTPUT  TAPE  JTP » 1012 ,OL2 ,E IGEN 
AVED*0. 

DO  231  I-1,N 

231  AVED-AVED+DIII) 

ATN»N 

AVED-AVED/ATN 

SATN»N-3 

IFISATN)232,232,233 

232  WRITE  OUTPUT  TAPE  JTP, 1018 ,AVED 
GO  TO  235 

233  STND»SQRTFIDL2/SATN) 

WRITE  OUTPUT  TAPE  JTP, 1019, AVE0,STND 
C  INTERROGATE  NEXT  CARD 

235  READ  INPUT  TAPE  ITP , 1002 ,N ,NO,L, I ZY 1 1) , I *1, 10 ) 

IFIN)240,240,10 

240  READ  INPUT  TAPE  ITP, 1013, N, LI, NO, L2 
1FIN)270,270,24S 
C  OUTPUT  TWO 


I 
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245  WRITE  OUTPUT  TAPE  JTP*1014 
GO  TO  260 

250  READ  INPUT  TAPE  I  TP. 1013.N»L1*NO.L2 
IF(NI2T0.270.260 
260  IF(Ln275.275.300 
275  IF(L2)280.280.340 
280  00  290  1-1.3 

XMKdl-VMKd.N) 

290  YMKdl=V^<Kd.NO) 

LLL-1 
GO  TO  350 

300  IF(L21310.310.330 
310  DO  320  1-1.3 

XMK(l)-UMKd.N) 

320  YMKd)-VMKd.NO) 

LLL-2 
GO  TO  350 
330  DO  335  1-1.3 

XMK( I )-UMK(  I .N) 

335  YMK(  n=UMIC(  I  .NO) 

LUU-3 
GO  TO  350 
340  DO  345  1=1.3 

XMKd)*VMX(  I.N) 

345  YMXd)-UMX(  I. NO) 

LLL=4 

350  DTPD=D0TPRD(XMX( 1) .XMK(2) .XMX(3) »YMX( 1) .YM<(2 ) . YMX ( 3 ) . B1 .B2 .B3 » 
ICOSA.COSB.COSC) 

WRITE  OUTPUT  TAPE  JTP.1006 
GO  T0(360.370.380.390) .LLL 
360  WRITE  OUTPUT  TAPE  JTP . 1016.N.NO.DTPD 
GO  TO  250 

370  WRITE  OUTPUT  TAPE  JTP . 1036.N.NO.DTPD 
GO  TO  250 

380  WRITE  OUTPUT  TAPE  JTP . 103S.N.N0.07PD 
GO  TO  250 

390  WRITE  OUTPUT  TAPE  JTP.1037.N.NO.OTPD 
GO  TO  250 
270  PRINT  1017 
CALL  EXIT 
1000  FORMATd2A6) 

1002  FORMAT(2I3.I6.rOA6) 

1003  F0RMAT(6E10.4) 

1005  F0RMAT(1H1.7X.12HPLANE  NUM3ERI 4, 1 10 ,6H  ATOMS .7X , 10A6// //12X .90HEQU 

2ATION  OF  PLANE  IS  Ml*X  +  M2»Y  +  M3*Z  =  D  WHERE  D  IS  THE  O'R  I 
3GIN  TO' PLANE  t)ISTANCE/lH0.20X.3HMl  =  E12.5.7H  M2-E12.5.7H  M3-E 

412. 5. 6H  D-E12.5////12X,32H':00RDINATES  OF  CENTROID  X  -E12.5 

5.5X.3HY  -E12.5.5X.3H2  •12.5////15X.72HATOM  NUMBER  WEIGHT 

6  0  DELTA  D  (DELTA  D )**2/ I IHO . 1 7X » I 3.6X . E10.4. 

75X.E12,5.5X.E12.5.4X.E12.5)) 

1006  FORMAT(IHO) 

1012  FORMAT(////15X.21HSUM  WT*(DELTA  D)**2  -E12. 5 .20X . 13HF IGEN  VALUE  -E 
112.5///) 

1013  FORMAT(4I3) 

1014  FORMAT! 1H1.29X.40HDIHEORAL  ANGLES  BETWEEN  PLANES  AND  L I NES////29X . 
134HANGLE  BETWEEN  PLANES  PI  AND  P2  =  A/1H0.2BX.33HANGLE  BETWEEN  LIN 
2ES  LI  AND  L2  »  B/1H0.28X.36HANGLE  BETWEEN  LINE  L  AND  PLANE  P  -  C) 

1016  FORMATdHO.28X.4HPl  »I3.i0H  P2  -I3.17H  COSINE(A)  -E12.5) 

1017  F0RMAT(1H1.9X.13HJ0B  FINISHED.//////////) 

1018  F0RMAT(1H0,14X.11HAVERAGE  0  -E12.5) 

1019  F0RMAT(1H0.14X.11HAVERAGE  D  -E12 .5 . 15X . 20HSTANDARO  DEVIATION  -E12. 
15) 

1020  FORMAT  dHl.35X.12A6) 

1021  F0RMAT(1H1.71H  A  SATISFACTORY  LEAST  SQUARES  FIT  CANNOT  BE  FOUN 
ID  FOR  SET  NUMBER  I3.38H.  EXAMINE  INPUT  PARAMETERS  FOR  ERROR.) 
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1U25  F0RMAT(1H1.8X»UHLINE  NUMBER  14 . 1 10 ,6H  ATOMS .7X , 10A6//// 12X ,40HTHE 
IPARAMETRIC  EQUATIONS  OF  THE  LINE  ARE//37X.4HX  *  E12.5.5H  +  E12.5 

2f4H  *  T//37X.4HY  =  E12.5t5H  +  E12.5.4H  ^  T//37X.4HZ  »  E12.5f5H 

3+  E12.5*4H  *  T////12Xt46HTHE  SQUARE  DEVIATIONS  FROM'  THE  LINE*  D** 

42*  ARE//  22X»11HAT0M  NUMBER  •  1 IX tSHWE IGHT . 14X «4H0**2/ ( 1H0*24X • I  3* 1 
54X*E12.5>6X*E12.S) ) 

1032  F0RMAT(////12X.21HSTAN0AR0  DEVIATION  »  E12.5) 

1035  FORMAT! 1H0*28X.4H).1  =I3»10H  L2  =I3»17H  COSINE(B)  «E12.5) 

1036  F0RMAT(lH0i28X.4H  L  -I3.10H  P  ai3,i7ri  SINE(C)  »E12.5) 

1037  FORMAT! 1H0.28X.4H  P  =I3,10H  L  =I3.17h  SINE!C)  =E12.5) 

END 


Point  to  Peak  Distance  Calculation 


liom _ DETBRMlMtB  PISTAMCES  OF  ATOMS  PROS  AMY  POINT  IN  UNIT  CELL 


DIMENSION  NAMeii21.Ht5Q0I.Xt50QI.YtB00>. _ 

lZ(SOOItNSll200)«NS2(200t«SV(200}tNl(SOOOI.N2(5000). 
2JNT(lIl.NINT(100.Il>.AVE(10QI.EONIIOO>.NTBtl001.Di5000>.TTIII>. 
3TVI1I)  .DISdDtOISAdOO.llI 

CQSPf(X)«CO$fia«IAWX/HQ«l _ 

I  TP-5 

READ  INPUT  TAPE  ITP*  1000.  (NANEd  1 1  * Il-l*  12) 

-R£AP  INPUT  TAPE  lTP.12QQiA.B.C»Ai..BEtftA _ 

READ  INPUT  TAPE  ITP.llOO.BNX 

READ  INPUT  TAPE  ITP.  1 lOO.XC.VC.ZC _ 

BNX-BNX*«2 


_ lifl _ _ 

100  I-dl 

REAP  INPUT  TAPE  ITP.llOO.Xd  >  .Yd  >  .!(  d _ 

105  Ndl-I-1 

_ IP(ABSP(XtHl»ABSP(Y(HI»A6SP(Zd)llllQ.U0.10Q 

110  JJ«0 


_ IlHl _ 

115  READ  INPUT  TAPE  lTP,1200.Tl.T2.T3tUl.U2.U3 

IP(ABSF(Tl)»ABSF(T2l4ABSP<T31»A6SPLUl»»A»SFtU2)4AB5F(U3) >140.140. 
1  120 

120  DO  130  J-l.l _ 

I1-1141 

N<II)-I1-1 _ 

Xdl>-Tl«Ul*XIJ) 

Yill)-T2»U2«Y1J> _ 

130  Zdl>«T34U3«ZIJ) 

CO  TO  115 _ 

140  CAL>COSOF(AL) 

CtE-COSOF(BE) _ 

CCA-COSOF(GA) 

150  00150  J-ltll _ 

157  V1-XU)-XC 

V2-YiJ)-YC _ 

V3-ZIJt-t£ 

VEC-ABSF<00TPRDiVl.V2,V3.Vl.V2,V3.A,B, _ 

IC.CAL.CB^.CCAIt 

_ IF(VEC-BNX)1B0.1B0,100 _  _ 

IBO  kK>KK4l 

N1IXK)-NU) 
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190  CONTINUE 

mEITE  Output  Tam  JT».»6A.(NANEU).K-l,12t 

00  200  J>ltKK 

Ll-*tl(Jl/l - 

M1-N1IJ)-L1«1*1 

200  MITE  OUTPUT  TAEC  JTR,200fl.Ll»Nl,b(Jt - 

READ  INPUT  TAPE  ITP.  llOOtXC»VC«ZC 

- lF(AISP(XCI«i¥SfTVCUAOSP(2Cn21«0.21«0«210 

210  KK-0 

- 60  TO  ISO - 

2190  CALL  EXIT 

1000  PORNATtllAAl _ 

1100  F0RMATI3E10.4) 

120Q  F0KHAT<6E10«9I _ 

2000  F0RNAT<lH0tlTX.12tlH-.I2t4X*P0.5l 

3300  FBRNATtlH1.10X.12AAI _ 

ENOdtliO.O.OtO.  1.0. 0.0.0.0.0.0.01 


Form  Factors  for  Busing  Least  Squares 

CFORMBUS  FORM  FACTORS  FOR  BUSING  LEAST-SO.  REFIN.  PROGRAM 
DIMENSION  S(33).F(32) .NC ( 25 ) .BLANK! 12 ) 

ITP«5 
JTP-6 
KT»10 
REWIND  KT 

READ  INPUT  TAPE  ITP. 250,  ( BLANK! I ) . I »1 ♦ 12 1 

NC! 1)«1 

NCI9I=2 

NC! 17)=3 

NC!25)=4 

5  READ  INPUT  TAPE  ITP. 

1  100.A.A1.B.61.C.ATOM 

IF!AI30.30.10 
10  S!l)*1.55 
DO  20  1-1,32 

F!  n=A*EXPF!-Al*S!  I  )**2)+B*EXPF!-Bl*S!  n**2)+C 
20  S! I+ll-S! I )-.05 

WRITE  OUTPUT  TAPE  JTP.  200.ATOM 
WRITE  OUTPUT  TAPE  JTP,  210 
WRITE  OUTPUT  TAPE  JTP.  220 
WRITE  OUTPUT  TAPE  JTP.  210 
WRITE  OUTPUT  TAPE  JTP. 

1  230.!S!  I  I.F!  I )  .S !  1416  >  .F !  I-MS ) » I  - 1 . 16 ) 

WRITE  OUTPUT  TAPEKT .250. !BLANK! 1) . I >1 . 12  I 

WRITE  OUTPUT  TAPEKT,240» !F! 1) ,F! I+l I ,F! 1+2) ,F! 1+3 ) ,F! 1+4) » 

IF! I+5),F! 1+6) ,F! 1+7 ) » ATOM.NC ! I ) » I *1 .25. 8 ) 

60  TO  5 

30  END  FILE  KT 
REWIND  KT 
CALL  EXIT 

100  FORMAT!5F8.5,A6) 

200  F0RMAT!1H1,9X.A6) 

210  FORMAT! IHO) 

220  FORMAT!1HO,9X.7HSIN! )/L.13X.llHF0RM  FACTOR, 29X,7HSIN! )/L,13X.llHF0 
IRM  FACTOR) 

230  FORMAT! 1H0.9X.F5.2.15X.F10.6.30X.F5.2.15X.F10.6) 

240  F0RMAT!7F9.4,F8.4,1HF,A6.I2) 

250  F0RMAT!12A6) 

END 
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Variance-Covariance  and  Atomic  Parameter  Input 
for  Busing  Function  and  Error 

CBUSVAR  VARIANCE  MATRIX  INPUT  FOR  BUSING  ERROR  FUNCTION  PROGRAM 
DIMENSION  A(18528)tBll92) 

ITP»5 

JTP-6 

I5«9 

REWINDI5 

READ  INPUT  TAPE  ITP. 

1  lOOiNiSF 
READ  INPUT  TAPE  ITP, 

1  200,A1,A2,A3 

DO  10  I»1,N,3 
READ  INPUT  TAPE  ITP, 

1  200,B(II,B(I-)-l),B(I-f2) 

5  B(  n«{B{  n/All**2 
BII+1)»(B(I+1)/A2)**2 
10  B(I+2|»(B(I+2)/A3)**2 
SF»SF**2 
M«0 

DO  40  I«1,N 
DO  40  J>1,N 
IF( I-JI20,30,40 
20  M»M+1 
A(M)«0,0 
GO  TO  40 
30  M=M+1 
A(M)«B(n 

40  CONTINUE 
MN*M/e 

IF(M-8*MN)45,45,41 

41  L«M+1 
M»MM*a+8 

DO  41  I*L,M 
43  A( I  )«0.0 

45  IFISENSESWITCH  1)46,47 

46  MM»1 
PUNCH  600 

;  PUNCH  100,N,SF 

•  GO  TO  48 

!  47  MM»2 

I  WRITE  OUTPUT  TAPE  1 5,100, N,SF 

I  48  K=0 

I  DO  55  I>l,Mt8 

I  K«K+1 

I  GO  T0<51,50)»MM 

I  50  WRITE  OUTPUT  TAPEI 5,300 , A ( I » ,A( I+l ) ,A( I4-2 ) ,A ( I +3 ) ,A ( 1+4 ) , A( 1+5 ) , 

I  1A(I+6),A(I+7),K 

I  GO  TO  55 

P  51  PUNCH  700,A(I),A(I+l),A(I+2),A(I+3),A(I+4),A(I+5)»A(I+6),A(I+7) 

f  55  CONTINUE 

I  READ  INPUT  TAPE  ITP, 

I  1  200,(B(ntB(I+l),B(I+2),I>ltN,3) 

f  IJ«N/8 

!F(N-IJ»8)80*80,60 
60  L-N+1 
N»8*IJ+8 
DO  70  I-L,N 
70  BII)-0*0 

80  GO  T0(81,82)*MM 

81  PUNCH  600 
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82  K-0 

DO  95  I-l»N>e 
K-K+1 

GO  TO (92 *90) (MM 

90  WRITE  OUTPUT  TAPEI 5 (SOOtB ( 1 1  tB ( I-^l )  »B(  1  +  2  I  (3 (  I ■t'S )  »B ( I')-4 )  tB (  U 5  )  t 
lB(I+6)«B(I+7)tK 
GO  TO  95 

92  PUNCH  800tB(I)«B(I+l)tB( l+2)«B<I+3)«B(l+4)«B(I+5)*B(l+6)tB(I+7) 

95  CONTINUE 
END  FILEI5 
REWIND15 
PRINT  400 
CALL  EXIT 

100  FORMAT(I9iE10i5) 

200  F0RMAT(3E10.5) 

300  F0RMAT(8E9.4«I4<3H  SO) 

400  FORMAT) 58H  JOB  FINISHED.  SAVE  TPE  A5  FOR  OFF-LINE  CARD  PUNCHIN 
IG. ///////////) 

500  F0RMAT(8F9-.6tI4t2H  P) 

600  FORMAT(72X) 

700  FORMAT(8E9«4) 

800  FORMAT(8F9.6) 

END 


Function  ARCSIN(X) 

FUNCTION  ARCSIN(X) 

10  A«+l«5707963 
20  B»-0«21459880 
30  C-+0. 08897899 
40  0»-0. 05017430 
50  E«+0.03089188 
60  F»-0. 01708813 
70  Ga+0. 00667009 
80  H—0. 00126249 
100  YY-X 

110  IF(YY)  120.210.120 
120  Y-ABSF(YY) 

130  IF(Y-l.O)  140.230.280 
140  2*-S<3RTF(1.0-Y) 

150  2»A+Z*(A+Y«(B+Y*(C+Y*(0+Y*(E+Y*(F+Y*(G+HY) )>)))) 

160  IF(YY)  170.210.190 

170  ARCSIN*  -Z 

180  RETURN 

190  ARCSINi>Z 

200  RETURN 

210  ARCSlN-0.0 

220  RETURN 

230  IF(YY)  260.210.240 

240  ARCSIN-A 

250  RETURN 

260  ARCSIN--A 

270  RETURN 

280  CALL  EN0J06 

290  END 


fflf 
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Subroutine  RECIP 


SUBROUTINE  RECI P ( AR iBRi CR i ALRtBER (GAR » AA«BB iCC tCOSALtCOSBE iCOSGA ) 

ALG»ALK*3.14159/lS0i0 

B£G«a£R*3tl4159/100.0 

GAG»GAR»3,14159/.18C.C 

C05AR=CCSF(ALG) 

CCSuR»COSF{B£G) 

C0SGR«C05F(GAQ) 

VR»A«*E>R*CR*SORTF( 1.0-COSAR*COSAR-COSBH*COSBR-COSGR*COSOR+2i*COSAR 
1*COS3R»COSGR) 

SINAR«5INF(ALG) 

SINSR»SINF(3EG) 

SINGR>S1NF(GAG) 

AA«3R*CR*SINAR/VR 

bB»AR*CR*SINBR/VR 

CC«AR»dR*SINGR/VR 

COSmL> ( COSOR«COSGR-COSak ) / ( S 1 NbR*SI NGR ) 
COSu£«(COSAR*COSGR-COSBK)/(SIi'iAR*SINGK) 

COSGA- ( COSAR*COSBR-COSGR ) / (S IN«R*S1NBR ) 

RETURN 

ENU 


Function  DOTPROD 

FUNCTION  OOTPRO(UtV«WiX«YiZiBltB2tB3>COSA«COSo>COSC) 
00TPR0«U*X*B1**2+V*Y*B2**2+W*2*33**2+(V*2+Y*W)*B2*B3*C0SA+ 
l(U*2+W*X)*Bl*63*C&So+(U*Y+V*XI*ol*e2*COSC 
RETURN 
END 


Function  SCAFAC 

FUNCTION  SCAFAC(AtAl«BtBliC»SI 

SCAFAC«A*EXPF(-Al*( S«*ZI l+B^EXPF (-B1»<S*»2 >  >+C 

RETURN 

END 


Subroutine  MTXMUL 


SUBROUTINE  MTXMULILtMtNf AfBtC) 

DIMENSION  A(3»3liB(3i3)iC(3f3l 
t  DO  10  I»lfL 

DO  10  J«1«N 

»  10  C(ltJI«0> 

!  00  20  I>1»L 

t  DO  20  J>1«N 

DO  20  K^lfM 

20  C(liJ>>C(I>J)-^A(I«K)«B(XiJ) 
s  RETURN 

end 

Function  INTG(A) 


FUNCTION  INTGIA) 

IF(A)20»10»20 
10  1NT6-0 
GO  TO  30 

20  lNTG-XSI6NFIXINTF(ABSFIA)«.l)*XlNTFtAI) 
30  RETURN 
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